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ABSTRACT 


Presently, little is known about the intensity of the general circulation 
of the atmosphere. If an accurate measure of the generation of available po- 
tential energy (APE) could be made then this intensity would be known. Using 
Lorenz' (1955) definitions, the volume integration of the covariance between 
the deviations of the diabatic processes and temperature from their isobaric 
area mean values is used to calculate this generation. 

A seven-layer quasi-nondivergent model is used to compute the net effect 
of the diabatic processes within each layer. The data are the daily objective 
height analyses for the Northern Hemisphere from the National Meteorological 
Center for March 196%. The surface area contained within the analyses extends 
meemeanout CON to the pole. 

Utilizing the simplified vorticity and the full thermodynamic equations, 
the lower boundary processes and the accumulated vertical changes in the total 
derivative (horizontal) of the geostrophic potential vorticity are used to cal- 
culate the diabatic effects. Effects of variable terrain height and surface 
friction are included. Static stability is assumed to be pressure dependent. 

Analysis errors apparently cause erroneous values to be calculated south 
of 40N. Also, the values above 40 cb are unusable. This is due to the failure 
of the quasi-nondivergent theory in this region. Therefore, only the results 
in 18% of the earth's atmosphere are available for analysis. 

Monthly mean averages of the diabatic processes in the layer from 100 to 


40 cb indicate that the transition of the surface cover to permanent snow-cover 


Ma 





is felt throughout the troposphere. Extensive heating areas at the 40-cb level 
imply that average high cloudiness is required in order to trap the long wave 
emission. 

Using the daily values of the diabatic effects, daily computations for the 
generation of zonal and eddy APE are made. In the troposphere in March 1963, 


-4 ~2 = 
the average generation of zonal APE was 45 x 10 kj m sec’ and eddy APE was 


destroyed at the rate of - 22 x nO. kj m sec, Therefore, on thevavcrager 
heating (cooling) occurred in the low (high) latitudes and the warm (cold) 
eddies were cooled (heated) for this month. If the remaining half of the hemi- 
sphere is included in future calculations, most likely the zonal generation 
value will be larger and the destruction of eddy APE will be smaller. 

For the region examined, most of the generation and destruction occurs in 
the lowest 30 cb of the atmosphere. There is a near balance between the zonal 
and eddy components in the lowest 15 cb which suggests that most of the net 
seneration of APE occurs above the boundary layer. 

Harmonic analysis reveals that waves 2 through 8 are responsible for the 
majority of the destruction of eddy APE. Waves 2 and 4 are the dominant modes 


which suggests that the land-sea distribution and thermal contrasts are quite 


important during March 1963. 
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CHAPTER 1 


ENERGETICS AND AVAILABLE POTENTIAL ENERGY 


fee ENTRODUCTION 

The concept of energy describes the ability of an object or a system to do 
work, i.e., to exert a force through a distance. Studies of the energetics of 
the earth's atmosphere are intended to show how the atmosphere works by deter- 
mining the energy flows among the various energetical forms. Energy forms 
Povemeare significant within the atmosphere consist of three distinet Ubypes: 
(1) kinetic energy, (2) potential energy, and (4%) internal energy. 

Kinetic energy is the measure of work that a moving object can perform in 
being slowed to a state of rest. An example is the energy the wind expends as 
it propels a sailboat. Obviously the faster a particle moves the more kinetic 
SHemey itt has. 

Potential energy defines the capability of doing work by virtue of an oOb- 
ject's position in an external force field. A mass at a greater distance from 
the center of the force field has more potential energy. Matter falling due 
tO the force of a gravitational field is an example of a body losing this form 
en energy. 

Internal energy may be defined as the sum of the electromagnetic potential 
Sweaty Or electrons, nuclei, and the molecules of matver and the kKinetre eneray 
iwe to the molecular motions. This energy is closely associated with the aver— 
age energy of the thermal motions of molecules, i.e., the temperature of an ob- 


ject. When two identical objects with different temperatures are placed in 





contact, the warmer mass transfers internal energy to the cooler body through 
molecular collisions until both objects have the same temperature or internal 
energy. 

Energy is important because all matter possesses it and because energy is 
emeetiserved property, i.6., it can be neither created nor destroyed buyeonl, 
transformed from one form to another or transferred to another system. Thus a 
system in isolation from the rest of the universe acts under the constraints 
of conservation of energy. Since the earth-atmosphere system shows no tend- 
encies for velocity and temperature changes over a long time it can be consid- 
ered to be an isolated system. This means that over a sufficiently long time 
the system must be in equilibrium with the environment. Since the system in- 
teracts with its environment only through radiational processes, insolational 
warming received from the sun must be balanced by a cooling loss through long 
wave emission. However these relationships are valid only over a long time. 
Substantial departures from equilibrium may exist for short periods. 

Energy studies based on atmospheric observations must employ the hydro- 
static assumption. This is true because the assumption is used in calculating 
these atmospheric data. The hydrostatic assumption states that the upward- 
directed component of the pressure gradient force is balanced by the downward- 
acting acceleration of gravity at all times throughout the atmosphere. A di- 
rect consequence of this assumption is that the vertical velocity of the system 
is not a free variable but must be determined from the other state parameters. 
Since the resulting vertical velocities required to maintain the balance are 


small the flow within the atmosphere can be considered to be basically horizontal. 





Within an hydrostatic atmosphere, Haurwitz (1941) has shown that the po- 
wemial and internal energy of a column are proportional. That 1s sane tuereoe— 
of one requires an increase in the other. Therefore it is convenient to com- 
bine potential and internal energy into one form. Total potential energy is 
the name given to the sum of these components in an hydrostatic atmosphere. 
Wiin-Nielsen (1968) shows that the energy flow for an atmosphere exhibiting 
this equilibrium can be displayed as in the upper diagram of Fig. 1. Here G(I) 
represents the generation of internal energy; c(I,K,,) signifies the conversion 
of internal energy into the kinetic energy of the horizontal flow; C(I,P) por- 
trays the relationship between internal and potential energy; D(K,) denotes the 
dissipation by friction of the kinetic energy of the horizontal flow. The direc- 
tion of the arrows is determined by noting that over a sufficient time period 
the generation and conversion terms must off-set the removal of kinetic energy 
Since the winds do not ultimately cease to blow. 

Another useful approximation available to the investigator is the quasi- 
geostrophic theory. Charney (1948) provided a system by which the magnitudes 
of the terms in the hydrodynamical equations could be compared through the use 
of nondimensional analysis. He formulated the results for the synoptic weather 
scale while Burger (1958) extended this concept to the planetary waves. The 
approximation consists of keeping only the terms of largest comparable magnitude. 
It gets its name because the geostrophic wind relationship is the resultant 
when the horizontal accelerations are neglected in the first two equations of 
motion. Use of the theory is made when the vorticity equation is formed by 


SPress=dliterentiatins the first two equations Of motion and subtraccine etn 








Fig. 1. Flow diagrams of energy for an atmosphere in hydrostatic equilibrium. 








geostrophic wind is not used in the divergence term in the vorticity equation 
and hence the "quasi" portion of the name. Results of numerical forecasts us- 
img the theory led to the problem of spurious anticyclogenesis. This is due 
to the divergent nature of the geostrophic wind where southerly flow is con- 
vergent and northerly flow is divergent. This pattern does not resemble the 
observed picture presented by observations. A stream wind obtained from a 
stream function was introduced in place of the geostrophic wind to correct this 
@emect. since the calculated stream wind is a horizontally nondivergent wind, 
the assumption will be referred to as the quasi-nondivergent thecry in this 
paper. Further, the H subscript on the symbols involving the horizontal wind 
meeroeped in the remainder of the text. 

Lorenz (1955) introduced the concept of available potential energy (APE). 
Heederined it as that energy which would be available for conversion to kineéwuic 
energy if the atmosphere were brought hydrostatically to a state of minimum 
total potential energy through isentropic readjustment. The minimum state is 
ieee berized by horizontal stratification of the potential temperature sur- 
faces. An energy flow picture using this definition is displayed as the middle 
diagram in Fig. 1. The letter "A" represents the portion of the sum of poten- 
tial and internal energy which can be transformed; G(A) represents the genera- 
tion of APE; C(A,K) is the conversion of APE to kinetic energy; D(K) denotes 
wiestrictional dissipation of kinetic energy. It is clear from this figure 
E@acein an hydrostatic, adiabatic, and frictionless atmosphere, the sum of 
avallable potential energy and kinetic energy is conserved. A decrease in one 


must be accompanied by an increase in the other. However APE is exactly zero 





when potential temperature surfaces are horizontally stratified. The kinetic 
energy of the geostrophic wind would likewise be zero in such a case because 
the pressure surfaces coincide with the potential temperature surfaces. Since 
wiew actual wind tends to be largely geostrophic in nature, large voliecec: 
kinetic energy must therefore correspond to large values of APE. Thus when 
there is a simultaneous increase of both energy forms, diabatic processes must 
be responsible. Using Lorenz' (1955) sample calculations, the average amount 
of total potential energy which is available for conversion to kinetic energy 
is about one percent. Oort (1964) indicates that the average kinetic energy is 
smaller than average APE by a factor of three-fourths while Wiin-Nielsen (1967) 
shows the factor to be one-half. It appears that there is sufficient APE for 
conversion to kinetic energy. 

Wiin-Nielsen (1967, 1968) has pointed out that knowledge of the intensity, 
1.e€., the rate of energy transformation, of the atmospheric general circulation 
is limited. If thorough measurements of G(A), C(A,K), or D(K) were accessible 
then this intensity could be ascertained. However the assessment of any of 
these quantities is not available as a well-known value for the globe or either 
hemisphere. Generation of APE requires knowledge of the distribution of the 
diabatic heating and cooling. Conversion of APE to kinetic energy is dependent 
upon knowing the vertical velocity. Dissipation of kinetic energy is based on 
imleclonal effects both in the boundary layer and in the remainder ei eine sa aics. 
phere. 

Lacking adequate knowledge about the horizontal components of the frictional 


force, D(K) could be evaluated as a residual term in the manner of Kung (1966, 





1967). Wiin-Nielsen (1965) notes that such a procedure is difficult. He 
shows that the calculation involves tendencies toward a balance between nu- 
merically large quantities. Thus extreme accuracy in the observations of at- 
mospheric data is needed in order to avoid errors. It is doubtful that such 
accuracy exists even in the dense data collection networks. 

The choice of computing either G(A) or C(A,K) is one of selecting equally 
Miieeave Options. Vertical velocities needed for the conversion term arc ycic 
Peomeuro: the divergent part of the horizontal wind. Forecast models employing 
the quasi-nondivergent and adiabatic assumptions have provided weather predic- 
tions for a little over a decade and the results indicate that these approxima- 
tions yield useful prognoses over short time periods. Based on the resolvable 
scale of the models this implies that the effects of the divergent wind and the 
diabatic processes are small over these short time periods. If adequate infor- 
Matton concerning the diabatic processes were available then the calculated ver- 
tical velocity fields using this knowledge would contain more useful information. 
This paper continues the examination of diabatic processes and the generation 


of available potential energy. 


1.2. PREVIOUS CALCULATIONS 

Many calculations for G(A) have used the approximate expression derived by 
Lorenz (1955). He first shows that the average APE depends on the variance of 
pressure relative to the isentropic surfaces. This is the result of keeping 
only the first term of a power series expansion for APE. He then approximates 


this variance as the variance of temperature on an isobaric level. Using the 





isobaric form, G(A) is computed by: (1) determining the deviations from the 
isobaric area mean for both diabatic processes and temperature, and (2) cal- 
culating the volume integral of the covariance of these deviations. Dutton 

and Johnson (1967), in criticizing the results of investigators who used the 
simplified expression, argue that the approximation neglects the generation of 
of APE due to heating (cooling) on an isentropic surface at higher (lower) pres- 
sure. (Isentropic surfaces exhibit strong upward slopes in a poleward direc- 
tion within the troposphere.) They also mention that the approximation does 
not retain dependence on total diabatic heating. However, Wiin-Nielsen (1968, 
Equations 51 and 56) shows that the approximation is energetically consistent 
with the quasi-nondivergent theory and must be used in studies where the theory 
is invoked. 

Following what has become a common practice, the distribution of any param- 
eter over the globe may be separated into a zonal or mean value and an eddy or 
remaining value. The eddies constitute the difference after the mean has been 
subtracted from the total. Oort (1964) compiled the results of many investi- 
gators and developed an annual picture of the atmospheric cycle in terms of 
zonal and eddy quantities. He gives a value of 31 x er = IJ m= sec~* for the 
generation of zonal APE and - 8x 10°* kj m* sec ~ for the generation of eddy 
APE in the space domain. (Values in the space domain are determined by per- 
forming calculations on each set of a time series and then averaging the re- 
sults. Using the terminology of Lorenz (1967), the results in the space domain 
answer questions concerning the standing zonal and standing eddy circulations. 


= 


Multiplication by 10° converts the units of 10°* kj m= sec’* to watts m=.) 





Oort's cycle implies that on the average there is heating where it is warm and 
cooling where it is cold along a longitude circle and just the opposite along 
a latitude circle. The bottom diagram in Fig. 1 shows the energy flow using 
the separation of APE and kinetic energy into zonal and eddy components. Di- 
rections of the arrows are based on Oort's (1964) results for the space domain. 
Smagorinsky, Manabe, and Holloway (1965) reported on the results of a nu- 
merical model which did not contain the effects of mountains, land and sea 


contrast or a hydrological cycle. They obtained a generation of 3 x Lon kj 


4 =2 


“3 kj m= sec’* for eddy 


m= sec’* for zonal APE and a destruction of -9 x 107 
APE. Manabe and Smagorinsky (1967) commented on a model with a hydrological 
cycle where the lower boundary consists of a wet and flat surface without heat 


* kj m= sec’* for eddy APE 


capacity. They computed a generation of 5 x 10° 
apparently due to the large release of latent heat in the tropics. 

Dutton and Johnson (1967) applied the expression for the generation of 
zonal APE in isentropic coordinates to a somewhat arbitrary assignment of the 
components of Lettau's (1954) heat budget within a two-layer model to compute 
a generation of 56 x 10°* kj m= sec *. Using the results of Kung's (1966) 
Study of the dissipation of kinetic energy, they concluded that the generation 
of eddy APE was positive and equal to 8 x 10°* kj m= sec™*. Kung (1967) ex- 
panded his data to include 1200 GMT observations as well as OOOO GMT soundings 
over North America. A new computation revised the value of dissipation from 


-1 a 
- This of course rea 


64 x 10°* kj m= sec™* downward to 41 x 10°* kj m= sec 
verses the picture proposed by Dutton and Johnson (1967) since the generation 


of eddy APE becomes - 15 x 10°* kj mo seca= using Kung's newer estimate. 
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Wiin-Nielsen (1968) has commented on the results of Dutton and Johnson 
(1967) and Kung (1966, 1967). He points out that the method employed by Dutton 
and Johnson applies to the mixed space-time domain while Kung's results refer 
to the space domain. (The mixed domain gives results which apply to the sum of 
the time-mean and standing zonal motions and the sum of the transient and stand- 
ing eddy circulations.) Oort (196) displays the dangers involved in mixing 
such computations in that the generation of eddy APE has opposite signs in the 
two regimes. Wiin-Nielsen (1968) also mentions the limitations involved in 


using cross-sectional data along a particular meridian to create zonal averages. 


moe tot PURPOSE OF THIS STUDY 

Results of previous investigations of diabatic processes and the generation 
of APE, other than numerical simulation, are based on models with limited ver- 
tical resolution. These models use observational or climatological data and 
examples are Wiin-Nielsen and Brown (1960), Brown (1964), and Dutton and John- 
son (1967). These investigations must confront the problems associated with 
sparse data coverage and analysis errors. Although the numerical simulation 
models do not suffer these defects they do impose certain artificialities due 
to ignorance of the real behavior of the atmosphere. To date they have not in- 
cluded a realistic treatment of the lower boundary effects. 

The present observational system of upper-air soundings of heights and 
winds which assumes the hydrostatic balance will remain the primary source of 
atmospheric data for some time to come. Application of the simple quasi-non- 


divergent theory seems reasonable in light of the amount and geographic distri- 
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bution of data. Eventually numerical simulation will offer results based on a 
more realistic treatment of the atmosphere. Then the findings of simple theory 
based on observed data can be used to evaluate the calculations of simple theory 
based on the results of numerical simulation. 

Mae purpose of this inquiry is to determine the feasibility of using a 
quasi-nondivergent model utilizing greater vertical resolution than previously 
used to compute the diabatic processes and the generation of APE. Results of 
a more realistic treatment of the lower boundary condition are also examined. 
Chapter 2 contains the development of the equations used to calculate diabatic 
heating and cooling and the zonal and eddy components of G(A). The data and a 
m~wemen ror ensuring vertical consistency of the analyses are discussed in Chapter 
4. Results of the computations are presented in Chapter 4 and the conclusions 
and a summary are provided in Chapter 5. The Appendix contains the finite dif- 


i eenee forms for the equations used in the various calculations. 





CHAPTER 2 


DEVELOPMENT OF THE EQUATIONS 


2.1. THE BASIC EQUATION FOR DIABATIC PROCESSES 

Wiin-Nielsen and Brown (1960) used a two-level, quasi-nondivergent model 
to calculate diabatic heating and cooling from observational data. Using the 
same model, Brown (1964) extended the study to a larger data sample and in- 
troduced a more reasonable lower boundary condition involving the terrain and 
frictionally induced vertical velocities. However, frictional forces were neg- 
tected in the vorticity equation and as in the previous investigation static 
stability was made a constant for the layer 80-40 cb. This paper is a continua- 
tion of that work with the added refinements of greater resolution in the ver- 
tical and a more realistic treatment of the lower boundary. In addition, the 
static stability is allowed to vary with pressure. 

The vorticity and thermodynamic equations for the quasi-nondivergent model 


including the effects of diabatic processes and friction have the form 


OF OF 
oe ee es 
Sa dA Gene) aa (cs 
and 
38 =) _ (2) _ | BH 
a EVE = + ow = - CF (2.2) 
Respectively. 


> 
In the above equations, V is the nondivergent part of the horizontal veloc- 


e e . = C ze ia 
ity vector which is computed from the stream function ¥, i.e., V = kxV¥; 


We 
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io OU 


- ~ 1s the relative vorticity; f is the Coriolis parameter where Poalg 


~ ox Sy 
, ; dp 
its value at 45N; w is the vertical p-velocity, at? where p is pressure and t 


> 


: One on 
= is the acceleration of gravity; V = a kas a Jj is the 


is time; g = 9.3 m sec 
> > 
horizontal part of the gradient operator where i and j are unit vectors; 0 = gz 
d( £ne) 


where z is the height above mean sea level; o = a is the static stabil- 


ity where a is the specific volume and 9 is the potential temperature; R = 287 


ae ~1 =i 


kj ton ~ deg ~ is the gas constant for air; a = 00h ky inom > deg ~ is the 
specific heat at constant pressure; H is the time rate of diabatic heating or 
cooling per unit mass where a positive quantity means heating; in and a are 
mieex and y components of the frictional force respectively. 

In (2.1) the vertical advection of vorticity and the terms expressing the 
turning of the vortex tubes have been neglected. This is equivalent to neglect- 
ing the vertical advection of momentum in the first two equations of motion 
Gromewhich (2.1) is derived. In the first term on the right side of (2.1), 
replaces (f+t€). Wiin-Nielsen (1959a) has shown that these deletions and ap- 
proximations are required in the simplified vorticity equation in order to 
Satisfy and be consistent with certain properties which the complete vorticity 


equation exhibits. 


Solving (2.2) for w gives the result 


1 - 1 o& RH 
7 --|2 Gs vo) <j ee 


If (2.3) is differentiated with respect to pressure, the equation 


* (22 ae = a oe ata ic) a 
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is obtained. Substituting (2.4) into (2.1) produces 


OF oF 
ot. OS(s v2(é8 2). 8 B 2(4) (= - <4) 
pee vi G*t) +5, ee 7 il ey ae (255) 


Letting 


: o(f 
4 - C a ee dp oO =) ) (255) 


represent the geostrophic potential vorticity, (2.5) can be rewritten as 


, - 
E - Vee De 
dp \op PR ot ; wee) 


Integrating (2.7) with respect to pressure from p = Po Where Lhe Subserlppne 
indicates a lower boundary value, to some arbitrary pressure p = p, gives the 


basic equation for calculating diabatic processes, 


(226) 





The terms inside the brackets can be computed from observations on pressure 
levels. Thus the diabatic heating and cooling at any level can be calculated 


provided that He can be determined. 


2-2. CALCULATION OF FRICTIONAL EFFECTS 

Phillips (1963) explained that a necessary condition for geostrophic bal- 
ance to be assumed is that frictional forces are small compared to the Coriolis 
acceleration. This assumption means that important frictional effects must be 


limited to the boundary layer. The average pressure at the earth's surface in 
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the Northern Hemisphere is about 9/ cb in the winter for the surface area north 
of 20N. For convenience it will be assumed that the effective lower boundary 
of the atmosphere is 92.5 cb. Assuming also that the depth of the layer in 
which the frictional forces vanish at the top is 15 cb, then the 85-cb level is 
the middle of this layer. The effects of the frictional forces in the boundary 


layer are assigned to the 85-cb level calculations. 


Let 
OT 
Ro SS (2.9) 
and 
or 
F = =o 2 ko 
: SS ( ) 


where tT and + are the surface stress components. Substituting (2.9) and 
x Ma 
(2.10) into the term representing the frictional effects in (2.8) leads to 
OF OF OT OT 
y see e y ; 


Soy) ee aie 


Since the stress vanishes at the top of the friction layer by assumption, using 


finite differencing (2.11) becomes 


5 OT oF =) 
ze : =) ~ ~Ap\ ox | oy - ae 
maere Ap = 15 cb. 


it the terrain Stress is further assumed to be proportional to the square 


of the surface geostrophic wind in the manner of Lettau (1959) and Cressman 
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(1960), these stress components can be approximated as 


= V = V d ; 
"x Sone os "y Dee O (2.15) 


Using (2.13), (2.12) can be written as 





Be. ole OT > BP, re) . 
Ap \ ox 51 : one CE OYo? - oy C52, = 


The equation of state, p = pRT, has been used in (2.14) to eliminate density 
which was itself assumed to be constant with respect to local space differentia- 
mem. In (2.1), P. is the pressure at the terrain height; om is the Coerfie 
cient of drag based on Cressman's (1960) work; i is the resultant geostrophic 
surface wind and Uo and ve are its components; Ty is the temperature at terrain 
height. The computation of pressure, temperature, and geostrophic surface wind 
euecume lower boundary is performed using techniques developed by Clarke and 
Lawniczak (1962) which utilize the terrain data of Berkofsky and Bertoni (1955). 
iitesprocedure is outlined in the Appendix. 

In (2.14) a further assumption is made that the local space variations of 


C - ageemefzligible. The final form for calculating friecticon is therefore 


D 
Ap\. ox Oy . ApRT ; 


where a is the surface relative vorticity. 





nf 


2.4. DIABATIC PROCESSES AT THE LOWER BOUNDARY 
; , OO 
Applying the hydrostatic equation, a = -qa, and the state equation to (2.2) 


yields 


Saas i WEE (2.16) 
or 


(2717) 





where oe is the static stability at the terrain height and oS is the verereall 
velocity in the lower boundary. Within quasi-nondivergent theory static stabil - 
ity can be at most a function of pressure and as such is modelled after Wiin- 


Nielsen (1959b) with the form 


{ 
{ 
i) 


: 


6 NY Se a) (2710) 


Ao) 


where o = 1 m= sec’* cb ~ and p = 100 cb. 


The omega in the lower boundary, , , is calculated following the results 


Lo 
of Cressman (1960) and Haltiner, Clarke, and Lawniczak (1963). Their results 


show that the vertical velocity induced by the terrain can be approximated as 
> 
Ww = V-Y (25s) 


and the frictionally induced omega can be computed from 


et 
Noa BPo“ploro (2.20) 
PE FRI. ; 
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where ©. is the terrain vertical velocity and D po is the friction vertical veloc- 


aty. Therefore 


= a ; : 
® 6 Oo + Dio (2.21) 


The integrand of (2.8) can be computed at the pressure levels for which 
data are available. Letting B equal the integrand and using the results of 
(2.17), the diabatic processes can be calculated at the required levels by us- 


ing the following summation: 


(2222) 





where BS is the result at the top of the ip-th layer. 


2-4. THE GENERATION OF AVAILABLE POTENTIAL ENERGY 

It has been shown by Lorenz (1955) that the computation of the generation 
ot APE, using pressure &s the vertical coordinate, depends on the ceverrtaice 
of the deviations of temperature and diabatic processes from their areal aver- 
ages. As shown by Wiin-Nielsen and Brown (1960), this generation or rate of 


energy conversion can be calculated from 


p=P 
- _R pod i( 2 
C= = I uF H (=) dSdp (2.23) 


C 0 
DP pee 


where G represents the generation of APE in the units kj sec™*; 5 1s the surtace 
area and dS = a“cospd) dep Where a is the radius of the earth, © is the lat- 
itude, and dA and dp are the increments of latitude and longitude respectively; 


the prime indicates that the value is the deviation from the area mean value. 
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The area mean value is calculated from 


2 
aa a ‘> en 
Cy) = = yF? J" ( ) cosparar (2.2h) 
8 “G1 “o 
Where the over-bar denotes the average of the parameter; 5S = 2na=(sin Mo-sin 3) 
is the area of the region and 1 and Qo are the limiting latitudes. 

Using techniques similar to those employed by Wiin-Nielsen, Brown, and 


Drake (1963), the generation term can be separated into wave-number domain such 


that 
G = 34. = ec (2225) 


where G. measures the contribution from the zonal or the north-south covariance 
and G. measures the eddy portion or the east-west covariance. Letting the 


Fourier series 


CO 


ao) 5 [Y (~) cosna + Z (p) sinnd] (2326) 


represent the diabatic deviation field and 


2) = y (9) fu Pt [y_,(¢) cosnA + z_(9) sinnh. | (2.27) 


represent the thickness deviation field, the contribution from the north-south 


covariance becomes 


Cale 
enRaX Oo l Qo 

= = d d ) ° 8 

G f = is YoY, cosdedp (2.28) 


i a 
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The calculation of the generation due to the eddies can be expressed as 


N 
1 /%e2 
+ : 
2, (Yy, + Z,2,] cospagdp (2.29) 


(2.30) 





and 


tRa<A poe a Jee 
eo y— FD [You a FZ ze il tcosae (259) 
nn nn J 


_— =] =i! 
* 1=1 O.P, J n 





where m determines the vertical resolution; k determines the latitudinal resolu- 
tion; N is the numerical limit on the harmonic components; and Ap is the lat- 


itude increment. Note that for a particular n, (2.31) becomes 





ntRa<A ap APs Jae 
G =, _ = ye 2, 2 COsd) ee ie) 
e,n 1=1 CDs gal im nen J 


Dividing (2.30) and (2.31) by § will give the generation in units of kj m~* 


sec ?, 





CHAPTER 4 


THE DATA AND THEIR USE 


4.1. SOURCE AND TYPE 

Data used in this investigation are the National Meteorological Center 
(NMC) routine objective height analyses which consist of the OOOO GMI and 1200 
GMI fields for the month of March 1963. They were obtained, along with the 
analyses from the Air Force Global Weather Central (GWC) for the same month, 
meomecne National Center for Atmospheric Research at Boulder, Colorado. For 
Pies parcicular month the NMC analyses at 100, 85, 70, 50, 30, 20, 15, and 10 
cb and the GWC set for 85, 70, 50, 40, 30, 20, and 10 cb were provided. 

The analyses consist of grid-point values from the NMC octagon containing 
1977 points. This array covers the Northern Hemisphere using a polar stereo- 
graphic projection true at 60N where the distance between grid points at this 
latitude is 381 km. The lowest latitude circle fully contained within the 
first two rows of boundary points is 2ON. Fig. 2 shows the surface area covered 
by this octagon. The right side is slightly shortened by the map limits. 

Wherever analyses were missing in the NMC set the analyses from the GWC 
set, if available, were used. Five fields out of 496 were absolutely missing 
and these five were formed by averaging the available analyses. Eight GWC 


fields were used to replace missing NMC fields. 


fees SIATIC STABILITY 


Static stability is quite important for the calculations involving (2.8). 
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As pointed out by Holl, Bibbo, and Clark (1963), the analysis of the atmospheric 
structure, in addition to satisfying the equation of state and the hydrostatic 
equation, must also be statically stable everywhere. Although deviations from 
static stability and hydrostatic equilibrium exist in the atmosphere it is 
doubtful that these deviations can be accurately represented by the data or 
should be contained in large-scale analyses. 

A 10-parameter version of the static stability model developed by Holl, 
Bibbo, and Clark (1964) in which temperature is linear in pressure to the BYE. 
power, was used to ensure the vertical consistency of the NMC analyses. Static 
stability is computed using the form 


-R/ 3 bee 


Pp ao 
CaS =RD p =a 


= (3.1) 
where o* has the units of velocity squared and is equal to op~. This model 
partitions the column from 100 to 10 cb into 8 layers. A one-to-one transform 
is developed between the heights at 10 mandatory pressure levels and the 100-cb 
height, the thickness for the layer from 100 to 50 cb and the static stability 
parameters for 8 layers. An inverse transform is available so that given one 
set the other can be generated. Fig. 3 presents in detail how the column is 
subdivided. Within the model, the heights at 100 and 50 cb are unaltered by 
the check because these levels receive the greatest amount of attention in the 
analysis scheme. The 100-cb analysis is closely related to the surface anal- 
ysis where the amount of surface data exceeds the available upper-air sound- 


ings by a factor of 8. Above 50 cb, there is less data generally available 


and the inakcuracies are larger. Also the 50-cb level incorporates as a first 
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Fig. 3. The separation of the atmospheric column from 100 to 10 cb into layers 
for use in the static stability check. 
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approximation to the final analysis a short-range forecast which is the result 
of special emphasis in the forecast model. However, the remaining levels are 
changed to make the computed static stability lie within a prescribed range. 
Table 1 shows the upper limit placed on the computed value for each layer. 


-2 


ehould this limit be exceeded, the limit minus 1000 m= sec” is substituted 


and the heights in the column are recomputed. Should the value be less than 


~2 


100 m* sec’~, LG 1s replaced with as value or 200 m= sec~* and the heights are 


likewise recomputed. 


TABLE 1 


UPPER LIMIT ON THE 
STATIC STABILITY VALUES 


Sigma 1 49000 
Sigma 2 25000 
Sigma 3 2 LOOO 
Sigma 4 23000 
Sigma 5 45000 
Sigma 6 35000 
Sigma 5 C00 
Sigma 8 45000 


3.3. RESULTS OF THE STATIC STABILITY CHECKS 
The 10-level model requires the heights at 100, 85, 70, 50, 40, 30, 25, 


20, 15, and 10 cb. Using the GWC 4O-cb analysis and linear interpolation to 
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form the 25-cb heights, a test was run on the 1200 GMI fields for 1 March. Re- 
sults are listed in Table 2 under the heading "Test A." The very high number 
of failures in the sigma 4 layer which contains the 40-cb level indicates that 
the GWC analysis is incompatible with the NMC fields. High failure rates in 
the sigma 3 and sigma 5 layers are the result of the dependence on the 40-cb 
heights. No further use was made of the GWC 4O-cb analysis and no tests were 


made on the GWC fields. 


TABLE 2 


NUMBER OF GRID-POINT FAILURES AS A RESULI OF 
STATIC STABILITY TESTS ON THE ANALYSES FOR 1200 GMT 1 March 1963 


Sigma Sigma Sigma Sigma Sigma Sigma Sigma Sigma 


Test ; = - rf : y - g Total 
A 70 119 707 ES foe -- -- a 1654 
B 76 oe) 56 16" 5172 a 12 O 14.08 
C 16 59 54 9 780 83 78 0 86) 
D 74 52 O 6 Si 7 56 O 23 


The next test involved using a Lagrangian three-point interpolation for- 
mula. The formula was applied to the heights at (0, 50, and 30 cb to obtain 
the 4O-cb height and to the heights at 30, 20, and 10 cb to compute the 25-cb 
heights. Table 2 list the results under "Test B." Improvement is clearly 
evident around 40 cb but the sigma 6 layer which contains the 25-cb level has 
unacceptable results. 

Next, a test using the Lagrangian scheme for the heights at 40 cb and a 


feicar fit at 25 cb was made. Test C” in Table 2 has the results and the 
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problem layer in this test is sigma 5 which contains the 50-cb level. In the 
calculation of sigma 5, the heights at 40, 40, and 25 cb are the dominant el- 
ements. Since the heights at hO and 25 cb are computed, it appears that the 
eheck is too sensitive. 

A tolerance limit was added to the check to lessen its sensitivity. The 
height value at 40 cb has the greatest effect on sigma 3 (excluding the layer 
containing the 40O-cb height) and the height at 25 cb likewise has the greatest 
mececy On sigma >. Prior to checking the column, a calculation of sigma 3 and 
Sigma 5 is made. Should either fail the check and a tolerance of + 15 meters 
on the appropriate height, either 4O or 25 cb, allows the check to pass, a + 
flag for that sigma is set. During the check of the colum should a sigma 
fail the test, the precomputed tolerance for the heights at 40 and 25 cb is 
added to the failed value. If the test then passes, the "fail flag" is removed. 
In addition, a consecutive check was programmed. If the sigma values in two 
adjoining layers should fail a special flag is set. The rest of the sigmas 
are computed without regard to the pass or fail criterion. Pressure heights 
are then recomputed incorporating the effect of the corrected value in the 
lowest sigma layer which failed. A complete stability check is then performed 
on the new height set. The results are listed in Table 2 under "Test D" and 


this became the final version of the static stability test. 


3.4. COMPUTATIONS OF H , B, H., AND G(A) 
p 
Fig. 4 shows the levels at which a and B can be calculated based on the 


available data. For convenience the lower boundary effects, represented by the 
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Fig. 4. The levels at which the integrand of (2.8), B, and the diabatic proc- 
esses, H, are calculated. 





eg 


G@iabatic term, a have been assigned to che 92.5-cb level seine ct teccs er 
Bie frictional forces are included in the calculation of B at the 85-cb level. 
Once H. fields have been computed the generation of APE can be determined. 
Based on the available data the vertical resolution parameter, m, in (2.30) 
Bid (2.51) takes on the value 7. The latitude-longitude increment is taken as 
2.5 deg and the region is set as that bounded by the latitude rings at 2ON and 


87.5N. Thus oi = 18.75N and oe = 88.75N while k takes on the value 28. 





CHAPTER 4 


RESULTS OF THE COMPUTATIONS 


iy olATIC STABILITY CHECK 

A static stability test discussed in Section 4.3 was performed on the 
Gata which consist of height analyses at 100, 85, 70, 50, 30, 20, 15, and 10 
cb. The time period covered runs from 1200 GMT 1 March 1963, to 1200 GMT 31 
March 1963. For these data the average number of grid points requiring a 
height correction in the vertical structure is 395. This number includes the 
effects of the combined NMC and GWC sets as well as the averaged fields which 
together involve four different analysis times. Height fields at 15 cb for 
OOOO GMT 6 and 7 March and 1200 GMT 19 March and the 100-cb heights for 0000 
GMI 7 and 15 March were formed by averaging the available analyses from the 
adjoining times. All of the NMC analyses are missing for OOOO GMI / March and 
the GWC set was substituted. GWC fields at 85 and 70 cb for OOOO GMT 15 March 
are also required. A significantly greater number of stability failures are 
encountered for these four periods as would be expected. Excluding their 
effect, the average number of grid-point failures is 369. Most of the stability 
problems occur in the sigma 1 and sigma 7 layers. An average of 160 failures 
in the sigma 1 layer and 104 in the sigma 7 layer are computed. Again, the 
total number of grid points is 1977 so about 19% of these points present some 
stability problem in the vertical structure. Now it is quite likely that these 
failures are random. No tests were made to determine otherwise. If in fact 


they are random, then time averaging procedures should remove their effect. 


40 





om 


However, since almost 20% of the vertical structure of the NMC analyses is 
Questionable relative to their static stability, numerical calculations based 
on these unchecked data would seem to be also questionable. 

Finally, Fig. 5 shows the comparison between the average of the computed 
Static stability for the month, indicated by the dashed line, and the analytic 
Bepresentation as indicated by the solid line. For the region of the Northern 
Hemisphere involved in this study the analytic form is quite obviously accept- 


able. 


4.2. MEAN DIABATIC PROCESSES FOR MARCH 1963 

Calculations of H and B were made using the data as illustrated in Fig. 
4, Computations were made using the finite difference equations given in the 
Appendix and the results are available at twelve-hourly intervals for the 
month. Fig. 6 depicts the zonally averaged values of the diabatic processes in 
the meridional plane as a function of latitude and pressure for March 1963. 
Values south of 3ON are not shown as the same unrealistic picture presented by 
Wiin-Nielsen and Brown (1960) resulted, i.e., the values rapidly decrease and 
become negative. This is probably due to a lateral boundary problem associated 
with the analyses which will be discussed shortly. The units are ens Kj sec”* 
ton”* and multiplication by 0.864 converts them to deg day™?. 

The stratospheric rates in the layer from 30 to 12.5 cb do not correlate 
with any known computations involving observational data or model atmospheres. 


For example, based on the numerical results of Manabe, Smagorinsky, and Strick- 


ler (1965) and Manabe and Hunt (1968), the values in this layer should be near 
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Fig. 5. A comparison of the monthly average of static stability based on data 


as indicated by the dashed line and the analytic representation as shown by 
the solid line for March 1963. 
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mO.6 x 10 ~ kj sec ton 

In contrast, the tropospheric detail in the layer from 100 to 30 cb appears 
comparable to previous results. Vernekar (1967) for example calculated the 
mean values for diabatic processes for January 1964. He calculated them using 
meecomputed vertical velocity forced only by the vertical variation of the hor- 
izontal eddy transport of zonal momentum. His results are displayed as Fig. 7 
and essentially show the same pattern from 100 to 30 cb. 

Excluding the results above 40 cb and south of 3ON, the area mean values 
for the layers from 100 to 30 cb are shown in Fig. 8. On the basis of the 
vertical nature of the zero line in Fig. 6, the latitude band centered at 50N 
was used to separate the heating and cooling areas. The computed values of H 
mee plotted as values for their respective layers. Overall the results indi- 
cate that there is little difference with pressure in the troposphere. Other- 


wese the maximum cooling occurs in the layer from 85 to 70 cb and maximum heat- 


fmiesin the layer from 50 to 30 cb. 


4.4. MONTHLY AVERAGE FOR THE VARIOUS LEVELS 

Within the troposphere and the lower stratosphere the important diabatic 
processes are the absorption of solar radiation, the emission of long wave 
radiation, the release of latent heat (provided that some condensate returns 
to the surface as precipitation) and the transfer of sensible heat across the 
earth-atmosphere interface. The method of computing the diabatic effect used 
here precludes the separation into contributions from the components. Some 


qualitative remarks can be made however. 
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Fig. 8. Area mean values for the diabatic processes as a function of pressure. 


The values are for the heating region and the cooling region as separated by 
the latitude circle at 51.25N. 
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Within the total depth of the atmosphere the cooling effect of the diver- 
gence of long wave radiation flux is the dominant element. In the middle lat- 
itudes and the tropics, the contributions from the net gain of latent heat re- 
lease, solar absorption and the boundary flux are important in the lower tropo- 
Senere, 1.e., the layer from 100 to 50 cb. Cooling effects of long wave radi- 
ation should dominate the layers above 50 cb in the high troposphere and low 
stratosphere and the polar atmosphere. However, in regions involving cloudi- 
ness, both radiational warming and cooling areas can be found as shown by John- 
son and Shen (1968) for example. 

Monthly averages for the lower boundary and the levels at 77.5, 60, 0, 
meee; >, ond 12.5 are presented in Fig. 9 through 15. The units are the same 
as in Fig. 6. Areas where average heating occurred during the month are cross- 
hatched. Values for maximum rates of heating and cooling are plotted in the 
appropriate regions. The latitude circle at 30N has been outlined to indicate 
the practical reas of this study. It also separates the hemisphere into halves 
relative to the surface area. 

Monthly averages of the diabatic processes due to the presence of a lower 
boundary are displayed in Fig. 9. Phillips (1963) has shown that a maximum 
heating rate consistent with the geostrophic assumption is of the order 0.1 kj 


-1 


sec ton~?, 


Within guasi-nondivergent theory there are also limitations on 
the frictional and terrain induced vertical velocities. Since the vertical 
velocities are included in the computation of diabatic processes in the lower 


boundary and are not used elsewhere it is sufficient to limit the magnitude of 


a -1 


the heating term. A limit was set at 0.7 kj sec” ton and an average of two 
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Fig. 11. Monthly average of the diabatic processes at 60 cb for March 1963, 


in the units 10°“ kj sec” ton”. 
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Fig. 12. Monthly average of the diabatic processes at 40 cb for March 1963, 


in the units 10°> kj sec”? ton”. 
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Fig. 14. Monthly average of the diabatic processes at 17.5 cb for March 
1963, in the units 10°* kj sec”* ton”. 
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Fig. 15. Monthly average of the diabatic processes at 12.5 cb for March 
1963, in the units 10 * kj sec.~ ton >. 
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erid-point values per analysis time failed the check and were reset to the limit. 

Patterns presented in Fig. 9 are based on calculations using (2.17). They 
are therefore the results of the averages of the local time change in and the 
horizontal advection of the temperature at terrain height and the vertical ve- 
locities induced by the surface. Effects of heat transport across the inter- 
face are probably wholly contained in this computation for the Lower boundary. 
The general picture of cooling over the land masses and warming over the ocean 
bodies seems to support this. Heating patterns in the western oceanic regions 
are well correlated with the warm Gulf Stream and the Kuroshio Current. 

The fact that this average for the lower boundary has a variable pressure 
base should not be forgotten. This base ranges from 104 cb over the oceans to 
4u7 ecb in the region of the Himalyas. Thus there are areas in which latent 
heat release and radiational warming due to low cloudiness are strong modifying 
influences in the lower boundary. 

Dickson and Posey (1967) have prepared snow-cover probability charts for 
the Northern Hemisphere. Their charts indicate the probability of the cover 
being 1 in. or greater in depth. Continental areas north of 50N have a factor 
of 0.5 for the month of March. These same areas north of 55N have a factor of 
1.0. Figure 6 shows that the change between the mean heating and cooling re- 
gions occurs at 53N in the lower boundary. Sellers (1965) has tabulated the 
distribution of various parameters in latitude bands of ten-degree width for 
the globe. Between 50N and 7ON the oceans cover about 36% of the surface area. 
In March, the oceans north of 7ON are either covered by ice or snow. There- 


fore the surface area north of 53N can be considered to be snow-covered for 
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this study. The vertical orientation of the zero line in Fig. 6 seemingly in- 
dicates that this effect is felt throughout the troposphere. 

In order to compute the diabatic processes at higher levels in the atmos- 
phere the effects of the variable lower boundary were assigned to the constant 
pressure level at 92.5 cb. Therefore the results in the other levels contain 
effects which are not really there. This is so because (2.22) is used to cal- 
culate the processes at all levels above the boundary layer. However, the num- 
ber of grid points at which the pressure at terrain height is generally less 
than 85 cb is small. Primary areas where it does occur are the Rocky Mountains, 
Bie sHimalyas, and the Greenland Plateau. 

Mean monthly 7O-cb height contours, taken from the June 196%, issue of 
Weatherwise, are depicted in Fig. 16. Long wave troughs are located along the 
mesvern coast of North America, the eastern Asiatic coast, in western Russia 
and the eastern North Atlantic Ocean. The middle latitudes are dominated by 
three waves from the western United States westward to western Europe while 
the flow pattern exhibits strong zonal characteristics eastward from the middle 
of the United States to the middle of the Atlantic Ocean. Diabatic fields at 
77.5 and 60 cb (Fig. 10 and 11) show that a good correlation exists between 
the extensions of the cooling regions into the lower latitudes and the trough 
positions. The effects of the minor trough near Hudson's Bay can be detected 
also. Andrews (1963) reports that record and near-record amounts of precip- 
itation fell in the Ohio and Tennessee Valleys during the month. Changes from 
cooling to heating in this region in the layers above the lower boundary are 


most likely due to the release of latent heat and radiational warming below 
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cloud bases. Andrews (1963) also reports that portions of the Southern Plains 
and lower Mississippi Valley suffered a precipitation deficiency. This area, 
particularly at the lower boundary and the 77.5-cb level, shows strong diabatic 
eooling. 

The remaining levels, Fig. le through 15, should contain mostly the effects 
of long wave radiation. London (1957) acknowledges that radiational warming 
beneath clouds and the base of the stratosphere does occur but probably is the 
exception rather than the rule. Extensive areas of diabatic heating in these 
figures imply that this is not true. However, an explanation requires the 
existence of multi-layer cloud decks and/or multiple tropopauses. But note 
that the large magnitudes of the heating and cooling rates, especially at the 
12.5-cb level, indicate that something is wrong with these calculations. This 
of course is the same conclusion that was reached from the results presented 
eric. 6. 

Results for the 40O-cb level appear to be acceptable. Besides the effect 
Of cloud decks, perhaps the effect of the release of latent heat above 50 cb 
is also required to explain these heating patterns. The answer cannot be re- 
solved by this investigation. However, the amount of water vapor generally 
thought to be present above 50 cb would not appear to be significant enough 
to offer a substantial contribution toward latent heat release. Thus the log- 
ical conclusion seems to be that the clouds are the ultimate cause of the heat- 
ing patterns. The presence of cirrus clouds would be quite effective in pro- 
ducing these warm areas as has been shown by Johnson and Shen (1968). At least 


the required cloud patterns needed to support the results in Fig. 12 would be 
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generally eastward of the trough positions shown in Fig. 16. Katayama (1966, 
1967a, 1967b) has carried out extensive calculations for the radiation budget 
of the Northern Hemisphere. Although many assumptions are made, particularly 
with respect to cloud properties, favorable comparisions of the total solar 
radiation at the earth's surface with other studies appear to justify his cal- 
culation procedure. He states that the "existence of high-level clouds has an 
intense effect on decreasing’ the net long wave radiation at 20 cb. Katayama 
(1967b) shows how low, middle, and high clouds can cause heating in various 
layers within the vertical structure through flux convergence of the long wave 
component. 

What appears to be the answer to the problems with the calculation at the 
higher levels can be seen in the patterns of Fig. 12 through 15. At the 40-cb 
level the analysis is beginning to breakup into small irregularities south of 
4ON. This feature is much more obvious at 25 cb. The results become so un- 
manageable at the other levels that no contouring is shown south of 40ON. Re- 
ferring to (2.8) and (2.22), the term representing the frictional effects is 
included in the computation of B at 85 cb only. Ignoring i which is of 
smaller order anyway, the remaining terms involved in computing the diabatic 
processes are the local time change in and the horizontal advection of potential 


vorticity. The computation at any level is the resultant of the vertical sum- 
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As a consequence of quasi-nondivergent theory, the terms representing the 
local time change and the advection must be of comparable order or time scale. 


Sample calculations for B show that this is true at 50 cb for example. But 
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for 50, 20, and 15 cb the advective UVerm is ereaver by sau slesse Ome tende moe 
magnitude in much of the boundary areas and at scattered interior grid points. 
Since the local time change of potential vorticity at 50 and 15 cb are of the 
Same order, the problem seeems to be with the horizontal advection term. 

An analysis of the computed potential vorticity using (2.6) indicates that 
the patterns are of type 1 as defined by Phillips (1963). This means that geo- 
strophic potential vorticity is characterized by "approximate uniformity." In 
other words, the values at any level do not exhibit a large range. This is 
true because of the quasi-horizontal nature of the isobaric surfaces. There- 
fore the wind used to advect the gradient of potential vorticity is suspect. 

The horizontal nondivergent wind used for advection is computed using the 


linear balance equation introduced by Shuman (1957). It has the form 


ae Yo vFP-wW (4.1) 





f fies 

Flisaesser (1968) found this to be the best approximation to use in calculating 
a nondivergent wind. The linearized form is the result of excluding the other 
terms in the divergence equation. Their exclusion is based on the fact that 
they are of order Rossby number cr smaller. The nondimensional Rossby number, 
Ro is the ratio of the characteristic wind or particle speed, U, to the prod- 
uct of the Coriolis parameter with a characteristic horizontal wave length for 
the motion, L. Within the atmosphere the patterns of interest are such that 


U < 100 m sec™* 


3 


L > 1000 km and f = 10°* see”*. Therefore Rs 1. The value 
1 is associated with minimum L and maximum U. For a characteristic atmospheric 


wind speed of 15 m sec-* and a characteristic wave-length of TO" m, RS = Ore: 
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One could argue perhaps that in the stratosphere the Rossby number is no Longer 
a small quantity and that discarded terms of order Bh must be retained in the 
equations. It is an observed fact that both the characteristic wind speed and 
wavelength increase with height in the troposphere and low stratosphere. JHow- 
ever, even allowing U to double and L to remain at’ 10° m, the Rossby number is 
only about 0.3. Thus the balance equation does not appear to be a bad assump- 
Sion in the layer from 30 to 10 cb. 

It should not be forgotten that much manipulation was involved in calcu- 
lating and correcting the static stability of this region. The sensitivity of 
the check implies that care must be exercised. Discontinuities of the lapse 
rate of temperature in a sigma layer due to the tropopause are probably the 
main cause of the problems. However, a new version of the check was programmed 
which requires only the data levels which are available. Height fields cor- 
rected by the new stability test and the unaltered original analyses for a 
sample period were used to recalculate the diabatic processes. The same fail- 
ures of the method were encountered above 30 cb with both sets. 

Significant errors in the analyses which would damage the calculations 
could be present. Mean heating for the layer from 100 to 30 cb in the merid- 
Senal plane as a function of latitude is shown in Fig. 17. The horizontal 
lines are the values based on an interpolation of the January and April tab- 
ulations of Davis (1963) for the 101 to 30-cb layer. Other comparisons can be 
made to the partial layer results of Wiin-Nielsen and Brown (1960) and Vernekar 
(1967) who includes the results of Brown (1964). These results essentially 


parallel the results presented here. Agreement in Fig. 1/7 seems reasonable 
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Fig. l{. Mean diabatic processes in the layer from 100 to 30 cb for March 
1963. Horizontal lines are the interpolated values from Davis (1963). 
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north of 5ON but south of 5ON there is a substantial difference. Holopainen 
(1967) questions the scheme of objective analysis employed by NMC on the basis 
of the results of Wiin-Nielsen, Brown, and Drake (1963, 1964). The reason is 
the deviation of their computations of the mean poleward flux of angular mo- 
mentum from other calculations he presents. Deviations are particularly pro- 
nounced south of 50N being larger by a factor of about 1.6 at 35N. 

The Fleet Numerical Weather Central (FNWC), Monterey, California, provided 
data for a sample period during the month under investigation. Analysis at 
FNWC is performed on a square grid containing 3969 points in which the equator 
is an inscribed circle. The NMC octagon is a subset of the FNWC prid. Fig. 
18 shows the difference when the NMC analysis is subtracted from the FNWC anal- 
ysis. The level is 50 cb, map time is 0000 GMT 5 March 1963, and the unit is 
meters. Boundary problems apparently imposed by the limits of the octagon are 
readily observable. in general, the NMC analysis is significantly higher on 
the boundary and the difference extends to 30N. If the differences in Fig. 18 
are characteristic of all NMC analyses, the poor results south of 30N may be 
explainable. Additionally, there are substantial differences in the top-half 
of the map. This is probably caused by a data sparsity in one analysis vs. 
the other. The dipole near the pole could be the result of a position dif- 
ference for a system. It should be also mentioned that there are differences 
in the techniques of numerical analysis as outlined by Cressman (1959) and 
Carstensen and Lawniczak (1966). 

Effects of analysis errors or uncertainties, whether caused by technique 


or lack of information, cannot be taken lightly. Zonal values for the diabatic 
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Fig. 18. Resultant differences when the NMC analysis is subtracted from the 
FNWC analysis for 50 cb. The map time is 0000 GMT 5 March 1963, and the unit 
is meters. 
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processes in the layer from 100 to 40 cb using three different analysis sets 
are shown in Fig. 19. The three sets are: (1) the original NMC analyses (NSS 
NMC), (2) the statically stable NMC fields (SS NMC), and (3) the fields pro- 
vided by FNWC (NSS FNWC). (The FNWC data were not checked for static stability. 
FNWC began using the vertical consistency scheme in June 1963, so the data used 
here probably do contain some inconsistencies.) There appears to be little 
merrcrence between the results using the different NMC sets. However, the re- 
sults of the FNWC data show striking differences. The outstanding deviations 
are in the polar region and the area from 20N to 60N. A requirement for addi- 
tional global data to eliminate analysis uncertainties could not be presented 
more clearly than in the preceding discussion. 

Pompolnted out previously, the calculation of diabatic processes at @lewer 
pressure levels in the atmosphere depends on the vertical summation of the B 
Bempubavion in layers at higher pressure. Thus once the results at a partic- 
ular level become erroneous they are propagated in the remaining levels at 
lower pressure. Phillips (1963) explains that for geostrophic motion of type 
im tne product of the square of the Rossby number with the Richardson number 
must be of order 1 for the characteristic scale of motion. The Rossby number 


has already been defined as 


R=—., (4.2) 


Wiin-Nielsen (1965, Eq. 2.10) shows that the Richardson number can be written 


as 
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Fig. 19. Mean values of the diabatic processes in the layer from 100 to 30 
ed using the three analysis sets as indicated. 





ari 


0 

RR. =—s 
i avy 
dp 


(4.3) 


This is true in a model using quasi-nondivergent theory where the static sta- 
bility is defined as stated in Chapter 2. Using (4.2) and (4.3), Phillips' 


criterion for geostrophic motion becomes 
UmapesONR cy 
f-12 oF = 
dp 


One of the critical terms in (4.4) is the vertical wind shear. Obviously 
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as the shear approaches zero from either side the relationship cannot be main- 
tained. This means that the quasi-nondivergent theory must fail in regions 
where a maximum or a minimum in the vertical wind profile exists or where U 
is independent of the pressure coordinate. The mean zonal wind as a function 
of latitude and pressure for March 1963, is plotted in Fig. 20 with the units 
msec -. Relationships between the maximum in the zonal wind and the diabatic 
patterns between 30 and 10 cb of Fig. 6 cannot be accidental. The greatest 
magnitude of diabatic heating in Fig. 6 is directly coupled with the jet stream 
in Fig. 20. Based on the discussion of the boundary problem around the sides 
of the octagon, it is reasonable to conclude that the magnitude of the zonal 
wind at and above 50 cb near the boundary is too large. But the shear char- 
acteristics associated with a jet still remain. 

The results displayed in Fig. 5 shows that the static stability parameter 
undergoes a change of at least two orders of magnitude in the layer from 100 
mo LO cb. In itself this is enough to invalidate (4.4) in the stratosphere. 


This range in the values for static stability is also included in the calcula- 
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tion of diabatic processes as can be seen in (2.22). Note that the vertical 
summation of B and the diabatic processes at the lower boundary are multiplied 
by a factor which has a range of one order of magnitude in the same layer. 

Phillips (1963) also shows that (4.4) is violated if the characteristic 
wavelength approaches the planetary scale. But there is a problem in assigning 
a value for the stability since his relationship is a function of density strat- 
ification. Increasing the static stability would tend to cancel the increasing 
wavelength in the denominator. Therefore, based on the preceding discussion, 
it can be expected that the quasi-nondivergent theory must fail where: (1) a 
level of maximum characteristic wind speed is encountered which causes the 
horizontal advection term to dominate the calculation, and (2) a region of sub- 
stantial increase in the value of the static stability parameter exists. Both 
of these conditions exist in the high troposphere so unacceptable results in 


this region using the method of this study must follow. 


4.4. THE GENERATION OF AVAILABLE POTENTIAL ENERGY 

Since the computations of diabatic processes above 40 cb seem unreasonable 
the values were not used in calculating the generation of APE. Also the values 
south of 3ON are questionable and they were not used. This means that the re- 
sults of the computation represent only 70% of the atmosphere which covers the 
surface area from 28.75N to 88.75N. Although (2.24) requires the volume in- 
tegration over the whole atmosphere, this is simply impossible at this time 
because of the lack of information over much of the globe. This inquiry there- 


fore measures the contribution from about 18% of the atmosphere's volume. As 
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for the exclusion of the stratosphere, Oort (1964a), Smagorinsky, Manabe, and 
Holloway (1965), Manabe and Smagorinsky (1967), Perry (1967), and Manabe and 
Hunt (1968) all show that most of the generation of APE occurs in the tropo- 
sphere. The reasonable results in the troposphere as shown in Fig. 6 justify 
making the generation calculations. 

In the application of (2.30) and (2.31), because of the limitations just 
mentioned, the vertical resolution parameter, m, is equal to 4 and the lat- 
mcoudinal band count, k, is equal to 24. The diabatic fields at 92.5, 77.5, 
60, and 40 cb are assumed to be representative for their respective layers in 
the same manner as in the calculation for H. Grid-point analyses for the di- 
abatic processes and the layer thicknesses were transformed into latitude bands 
of 2.5 deg width containing 256 points using bi-linear interpolation. These 
bands were then separated into Fourier components using harmonic analysis and 
only the first fifteen wave components and the mean were used in the calcula- 
tions of the generation of APE. 

Diabatic and thickness values are available at twelve-hourly intervals 
from OOOO GMI’ 2 March through OQOO GMT 31 March. Likewise, the computations 
for the generation of APE cover the same period with the same interval. Indi- 
vidual results for each set were then averaged over the period in order to 
obtain the means. 

4.5. THE GENERATION OF APE AS A FUNCTION OF LATITUDE FOR THE LAYER FROM 100 
TO 30 cb 
The monthly average of the generation of zonal APE as a function of lat- 


itude for the layer from 100 to 30 cb is plotted in Fig. 2l where the units 
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Fig. cl. Monthly average of the generation of zonal APE as a 


function of 
latitude for the layer from 100 to 30 eb for March 1963. 
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are 10°* kj m* sec * and the values are related to the areas of their respec- 


tive latitude bands. Each band's contribution involves an area weighting fac- 
tor which obviously has a greater value in the lower latitudes. For example, 

ene area for the band centered at 3ON is larger than the area for the band at 

SON by a factor of about 1.7. 

An obvious result shown in Fig. 2l is that the plotted values are mostly 
positive. This means that on the average the zonal processes are positively 
correlated. Diabatic heating therefore takes place in the warm low latitudes 
and diabatic cooling occurs in the colder polar latitudes. A small negative 
region between ON and 50N exists because the transition zones for the area 
deviations of the mean diabatic processes and the mean temperature field do 
not coincide. Based on the relationship between the values shown for the area 
from 30N to 60N, it would be easy to imagine a mirror image extending from 45N 
to about 7.5N. Considering the area factor this would mean a substantial ad- 
dition to the generation of zonal APE. 

Average monthly results for the generation of eddy APE are shown as Fig. 
22. Here an obvious result is that the values are all negative which means 
aac, on the average, warm eddies are cooled and the cold eddies are warmed 
everywhere. Maximum destruction of eddy APE occurs in the latitude band cen- 
tered at 65N while a secondary maximum exists at 50N. Brown's (196) results 
for April 1961 and 1962, also show peaks at 65N and 50N. He likewise computes 
minima in the destruction of eddy APE between 7ON and 75N as well as near 5/7N 
and 4}5N. A probable explanation for the maximum destructive zones in Fig. 22 


1s that they represent average frontal positions and therefore the areas in 
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Fig. 22. Monthly average of the generation of eddy APE as a function of 
latitude for the layer from 100 to 40 cb for March 1963. 
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which maximum baroclinic activity occurs. Weather patterns associated with the 
arctic front are related to the area between 60N and 7ON. Intense circulation 
features characterized this band throughout the month as indicated by the hem- 
ispheric surface analyses provided by FNWC. The most persistent feature was 
the low-pressure center near Greenland and Iceland. Long wave radiational and 
boundary cooling of the warm eddies would be quite effective in destroying 
eddy APE along this front. Air masses with temperatures lower than the freez- 
ing temperature would also contribute to the destruction if they moved over 
surfaces with temperatures closer to the freezing mark. Daily values indicate 
that destruction occurred 91% of the time in the rings centered at 62.5N and 
65N. Significant periods of generation were calculated for the remaining rings 
which helps to account for the relative maximum of destruction at 65N. Polar 
front activity is responsible for the patterns in Fig. 22 between 3ON and 5ON. 
Precipitation processes and the resultant release of latent heat in the warm 
eddies coupled with the trapping effects of clouds are most likely the reason 
for the variability shown. These processes would of course reduce the net 
destructive effect. 

Based on observations in the latitudes south of 3ON, it is doubtful that 
destruction or generation of eddy APE could be substantial in that region. 
The familiar wave patterns of the middle and high latitudes simply do not exist 
there. Effects on the scale of short waves may be important in the tropics 
but the lack of adequate observations plus the present resolvability of the 
numerical grid prevent any meaningful calculations to be made. 


Standard deviations for each latitude band are tabulated in Table 3 along 
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with the values plotted in Fig. 21 and 22. Variability about the mean on a 


day-to-day basis is of course indicated by the magnitude of the standard devi- 


TABLE 3 


LATITUDINAL VARIATION OF THE GENERATION OF ZONAL AND 

EDDY AVAILABLE POTENTIAL ENERGY IN THE LAYER 100-30 cb 

COMPUTED FOR THE RING CENTERED AT THE INDICATED LATI- 
TUDE IN UNITS 10°* kj m7 sec™* FOR MARCH 1963 


Standard Standard 
eas ae Deviation - Deviation 
30.0 62 .49 63 06 =lo063 22.72 
Bee: 46.67 MORE -20.90 235 54 
SAO. 31.88 29.77 -22 44 40.82 
B55 22.02 22 .66 es lo) 3h 83 
40.0 Sun 7.3 -2h 17 50255 
ees. =5).07 INS a19 556 Bowie 
S56) -7 .20 i.e =e, gs 
7.5 Sy ih 26.60 -18 .30 46.90 
OO) sy Ie 28 .12 -26.61 41.87 
pa D Slee 559 -20.49 BS) Oy 
DE. O 45.40 Hip oe =e) he: 40.35 
Sie Oe] 68.13 = ae 41.28 
60.0 60.20 fe Fie Ie pile ée. 
b2 55 Woe ae alO ney -h6 .42 47.12 
65.0 100.78 LO S525 59, oO .33 
Sis, oon oo ne =35.90 Dero 
HO 0 100-50 Gi 55 SI Se 42.06 
(eS Alina 123.86 =O) NS eer. lak 
Fase, aie 192 .28 = 1165 10 Fishes 
Fir aes. 116 .06 247 .39 See 66.12 
Cen i 72 261.12 - 6.69 54 .99 
S255 167 .86 266 44 OL 95 41.34 
apy e: te al 412.26 OL oe 41.06 
OeD 89.14 426 .73 =o 13 .66 


ation. Weather phenomena are quite changeable about the globe and must nec- 
essarily produce such variations on a daily basis as recorded in Table 3. 


Assuming normal distribution, generation of zonal APE occurred most of the 
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time between 30N and 4ON and from 55N to 72.5N. Between 4ON and 5ON the tran- 
sition zone is quite variable. North of 7ON the standard deviations indicate 
that substantial changes occur over the period examined. Changes in the weather 
within the small area of the polar cap certainly contribute to this variability. 
Because of the small area and the characteristic scale of weather patterns, 
most of the generation of APE is confined to the zonal component. Daily values 
for the generation of zonal APE reveal that north of 7ON significant periods 

of destruction are encountered throughout the month. The length of these pe- 
riods increases in the bands closer to the pole. Since centers of low pressure 
and high pressure exchange positions around and at the pole several times 
throughout the month, the wide range of values per area in this small region 

is not surprising. Now it might seem rational to fault the numerical analysis 
in this area and blame the large standard deviations on the lack of data. 

While it is true that significant data coverage probably terminates near 7)5N, 
the very nature of numerical analysis prevents the area north of that latitude 
mom exhibiting a random character. The first approximation for the analysis 
is generally a short-range forecast in which persistence of the existing pat- 
tern is most easily recognized. If the area were completely devoid of informa- 
tion then this persistence of the first approximation would control the final 
version of the analysis. Variability about the mean would be small as a re- 
Sult. However, the analysis about this region would have to be essentially 
unchanging in order for that to happen. Examination of the results at 87.5N 
shows that the average length of either a period of destruction or generation 


of zonal APE is about 54 hr. Therefore the conclusion is that the variability 
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is in fact real. 


4.6. GENERATION OF APE WITHIN PRESSURE LAYERS 
A comparison of the average generation of zonal and eddy APE for the month 
as a function of the pressure layers used in the calculation is shown in Fig. 


“1 and the area 


23 and 2 respectively. The units are 10°° kj m= sec™” cb 
covered extends from 28.75N to 88.75N. The computations indicate that in the 
troposphere the largest contributions to these generations are found in the 
lowest 430 cb of the atmosphere. The generation of zonal APE is practically a 
constant in this layer where the contribution slightly exceeds that of the 

next hO ch. Significant destruction of eddy APE obviously occurs in the lowest 
40 cb probably due to the long wave radiational cooling of warm air masses and 
the heating (cooling) of cold (warm) air masses by the sensible heat transfer 
across the earth-atmosphere interface. Destruction of eddy APE decreases with 
height ee ting that the effects of latent heat release and differential long 
wave radiational cooling are involved. Clouds in the warm air would lessen 
aeweocoling effect of the long wave component and as mentioned earlier, there 
even could be substantial warming in the layers beneath these clouds. Strong 
cooling of clear and cold air masses coupled with the release of latent heat 

in the warm sector would act to lower the net value for destruction. Finally, 
the layer from 50 to 30 cb does show a positive generation for the eddy com- 
ponent. 


-1 


Suomi and Shen (1963) calculated a value of 5.8 x 10°* kj m= sec™~ for 


the generation of eddy APE in the layer from 100 to 10 cb using Explorer VII 
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Fig. 23. Monthly average of the generation of zonal APE as a function of 
pressure layers for March 1963. 
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Fig. 24. Monthly average of the generation of eddy APE as a function of 
pressure layers for March 1963. 
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radiation data for a portion of the Northern Hemisphere. Corcoran and Horn 
(1965), using Tiros 2 radiation data and synoptic data from western Europe 
westward to Japan between 3ON and 5ON, computed a mean generation of 0.55 x 
10°* kj m= sec”™* in the layer from 100 to 50 cb. Johnson (1967) examined 198 
radiometersonde flights from thirteen stations in the eastern United States 
and the Caribbean. He computed an average destruction of eddy APE of - 0.2 x 
io * kj m= sec™* due to long wave radiation in the layer from 60 to 2% cb. 
Clearly there is a large variability in the calculations which suggest that 


they are strongly dependent on time and scale and perhaps method. The standard 


deviations listed in Table 4 for the contributions from the layers indicate 


TABLE 4 


LAYER VARIATION OF THE GENERATION OF ZONAL 
AND EDDY AVAILABLE POTENTIAL ENERGY 
FOR THE REGION FROM 28.75N TO 88.75N IN THE UNITS 


10°© kj m= sec”* cb~* FOR MARCH 1963 


Standard standard 

mee a Deviation an Deviation 
100 - 85 cb oO 62 .9h =e 40.94 
S5) Ss We els Faron 38.11 Bip iy le eh .76 
{O - 5O cb 64. .69 ot ot Uo SA 255) 
SOe=2 50 eb A TAO) erat 4% 44 Peat 


that the most variable processes relative to the means are involved with the 
eddy APE above 70 cb. An examination of the daily values reveals that the 
generation of eddy APE was negative all of the time in the layer from 100 to 
85 cb. Likewise, the same is true for the layer from 85 to 70 cb with the ex- 


ception of one analysis. Occasional periods of generation existed in the layer 
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from 70 to 50 cb while the longest such period lasted from 1200 GMT 28 March 
through OOOO GMT 30 March. Results in the layer from 50 to 40 cb show that 
during the period from 0000 GMT 10 March through 1200 GMT 23 March, significant 
generation of eddy APE occurred. A shorter period from 1200 GMT 27 March 
through OOOO GMT 30 March was also characterized by generation which parallels 


the results in the layer immediately below. Table 5 lists the contributions 


TABLE 5 


INDIVIDUAL WAVE CONTRIBUTIONS TO THE GENERATION OF 
EDDY AVAILABLE POTENTIAL ENERGY IN THE LAYER 
FROM 50 TO 30 cb IN THE UNITS 10°* kj m= sec™* 
FOR THE INDICATED TIMES 


1200 GMT 1200 GMT O000 GMT 1200 GMI 

7 March 12 March 26 March 29 March 
Mega r 

Generation 2027 sie; =leae 12 .46 

Wave 1 -1.9 Oza Cine 2 
Wave 2 -1.9 ee -4..9 O.8 
Wave 4 amit -0.0 clea Del 
Wave O56 eee 2 6 2009 
Wave 5 -0.2 Cay => .4 O23 
Wave 6 -2.h 1.8 -1.9 0.6 
Wave 7 See lie -O.h oat 
Wave 8 oi, Al Lee -O.4 Oral 
Wave 9 Ce 0.6 LS =0.2 
Wave 10 Zon =OnG =i, 52 One 
Wave ll malee.) -O.4 -Q.4 -O.1 
Wave le -Q.2 O.9 QO.2 O.O0 
Wave 13 One -O.1 -0.2 0.4 
Wave 14 =O -0.7 Ope -0.0 
Wave 15 -O.1 OFe O20 -O.1 


from the 15 waves in the layer from 50 to 40 cb for four selected times. 


were selected on the basis of having the largest magnitude within a period 


They 


2 


@aaracterized by either generation or destruction.j This) magnitudes lc ecamered 
on the first line in the table and is practically the same for the four periods. 
Underlining has been used in the table to mark the wave component with the 
largest contribution within each set. Apparently the processes involved in 
waves 2 through 4 are most important in the generation or destruction of eddy 
APE in the middle and upper troposphere. This of course would not be true if 
the remaining waves acted in an opposite sense and were of sufficient magnitude. 
But in the data displayed, the first nine waves which contain most of the con- 
tribution tend to act in the same sense. Why the waves act so from one time to 
another cannot be answered here. The required relationship between the heat- 
ing and cooling regions and the temperature patterns to produce either genera- 
tion or destruction is quite clear. A cloud analysis for this layer might 
prove useful in explaining the results. 

A closer look at Fig. 23 and 24 shows that there is a near balance in the 
lowest 15 cb between the generation of zonal APE and the destruction of eddy 
APH. Excluding the lower boundary therefore would not detract from the calcu- 
lation for the generation of total APE or the intensity of the general circu- 
lation. Obviously the intensity of the zonal and eddy components would be mis- 
leading if this layer were not included. Of course these results are for one 
month only and for a limited region. Extensive conclusions should not be re- 
cited. But based on these averages, it appears that the diabatic processes at 
work in the lowest 40 cb of the atmosphere are the dominant factors in genera- 


ting or destroying available potential energy. 
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4.7. GENERATION OF EDDY APE BY WAVE NUMBER 

Fourier coefficients calculated by harmonic analysis were used to deter- 
mine the importance of the individual waves in generating or destroying eddy 
APE, Figure 25 presents the average of the daily values for the layer from 
100 to 30 cb in the units of 10°> kj m™* sec™* for the area from 28.75N to 
88.75N. With the exception of wave 1, Fig. 25 resembles the January 1963, 
picture presented by Brown (1964). Whereas he indicates that wave 1 is quite 
destructive of eddy APE for that January, these results show that the contribu- 
tion from this wave is not substantial when averaged for the month of March. 
However, Table 6 which lists the monthly averages along with the standard de- 
viations, shows that wave 1 has large variability. Obviously it is the most 
important single wave when it comes to substantial generation of eddy APE. 
Analysis of the daily values for the layer shows that 54% of the time the con- 
tribution from wave 1 was positive. In contrast, waves 2 and 3 generated eddy 
APE only 22% of the time. 

The dominant modes involved in destroying eddy APE within the layer are 
waves 2 through 8. Although the shorter waves appear to have a small role, 
lack of significant data over much of the hemisphere prevents a meaningful 
analysis of the importance of these waves. Waves 2 and 3 account for about 
50% of the destruction of eddy APE and apparently demonstrate the importance 
of the land-sea distribution in the Northern Hemisphere in forcing preferred 
modes. At least for this month, the cooling of warm eddies and the heating 
of cold eddies due to the scale of the land and ocean contrasts appear to be 


important effects. Figure 26 through 29 contain the monthly averages of the 
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TABLE 6 


HARMONIC ANALYSIS OF THE GENERA- 
TION OF EDDY AVAILABLE POTENTIAL 
ENERGY IN THE LAYER 100 - 40 cb 
FOR THE REGION 28.75N to 88. 75N 
IN THE UNITS 10°* kj m= sec™* 
FOR MARCH 1963 


Wave iene standard 
Number Deviation 
1 -0.21 4.00 
2 -5 34 oPale: 

3 -4..77 Dea? 

y 25518 4.90 
5 =% 28 Boe 
6 22515 pot 
i «0.81 2.74 
8 -1.03 1.80 
9 -0.00 le 
10 ~0.5), el 
ai -0.22 0.93 
12 <6) 1015 0.79 
2 0.06 One 
1h =o705 OF65 


15 =), 0.49 


latitudinal and harmonic analyses for each pressure layer used. Several fea- 
tures stand out in these figures. Waves e2 and 3 are destructive in all the 
layers while wave 1 and the medium scale waves gradually shift to a generating 
mode with decreasing pressure. Next, the baroclinic zone associated with the 
arctic front is pronounced throughout the troposphere while the polar zone 
shows the modifying effects of latent heat release and cloudiness with de- 
creasing pressure. The magnitude of the generation of eddy APE in the layer 
from 50 to 40 cb centered at 45N explains why the net result for the layer is 


characterized by generation. 





76 


1O7 KJ M2 SEC" 





100-85 CB 
90 80 70 60 90 40 30 
LATITUDE 
l 2 3 4 9 6 7? 86 9 HOI wicwo aia tS 


1O> Ky MSEC"! 


lOO-85 CB 





lb. 2. S3.°4 3 86. [7 WSS IO ime Visa 15 
WAVE NUMBER 


Fig. 26. Monthly averages of the latitudinal variation and the individual 
wave contributions for the layer from 100 to 8 cb for March 1963. 
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Same as Fig. 26 but for the layer from 85 to 70 cb. 
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Same as Fig. 26 but for the layer from 70 to 50 ecb. 
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Fig. 29. Same as Fig. 26 but for the layer from 50 to 30 cb. 





80 


4.8. THE NET RESULTS OF THE GENERATION OF APE 

In a study of this nature, many hours could be spent in exhaustive examina- 
tion of the results. But the data limitations and the small sample size only 
justify the highlights presented in this chapter. As mentioned in Section 4.3., 
a sample period based on different analyses was made available by FNWC. Table 


{ lists the results of the generation of APE calculations using the three dif- 


TABLE 7 


COMPARISONS OF THE GENERATION 
OF AVAILABLE POTENTIAL ENERGY 
BASED ON DIFFERENT ANALYSIS 
SETS FOR 0000 GMT 5 MARCH 1963 
FOR THE REGION 18.75N to 88.75N 
IN THE UNITS 10°* kj m7™* sec™* 


Analyses Gz Ge 

SS NMC Ons -38 2h 
NSS NMC hug =910) (1: 
NSS FNWC Lely er =e 


Berenu data sets. Recall, as shown in Fig. 19, there were not very marked 
differences between the two NMC sets relative to the mean diabatic processes 
in the layer from 100 to 30 cb. But the generation computations show how sig- 
nificant the subtle differences are. Zonal generation of APE is more than 
doubled and slightly more destruction of eddy APE occurs with the statically 
stable data. The difference between the FNWC and the NMC results relative to 
the zonal generation is disturbing. For this reason, any detailed investiga- 


tion of the daily calculations at this time is unwise. It may be that the 
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average results have a bias but that remains to be determined by performing 
the calculations on different data sets. With these thoughts in mind, the 
monthly average of the daily values for the generation of zonal APE is 45 x 
10-* kj m7? sec™* with a standard deviation of 20 x 10° * kj m= sec *, and for 
eddy APE is - 22 x Ore kj m-= sec7* with a standard deviation of iby oe Os k3 


-2 -1 
m eC . 


4.9. CALCULATIONS FOR JANUARY 1969 

Additional calculations have been made on data for January 1969. Analyses 
at 100, 85, 70, 50, 40, 30, and 20 cb were provided by FNWC which consist of 
grid-point values at 3969 locations for each level. The data are complete, 
Pimedliy consistent, and cover the entire Northern Hemisphere. 

Input observational data used by the FNWC objective analysis scheme are 
collected via a computerized network operated by the U. S. Air Force. Collec- 
bion centers at High Wycombe, England, and Fuchu, Japan, forward weather data 
tO Tinker Air Force Base, Oklahoma City, Oklahoma. These data are then for- 
warded to domestic users. Under normal operating conditions, more than 4500 
surface observations are available within 4 hr after OOOO GMT and 1200 GMT and 
about 400 upper-air soundings are reported within 4 hr of these times. 

FNWC reruns the upper air analysis using all data received within 10 hr 
of observation time. These data also include the results of a tropical anal- 
Bets COllected by a computer link to Fleet Weather Central, Pearl Harbor, 
Hawaii, which consist of winds at 5 deg latitude and longitude intersections. 


The region covered by these winds extends from 10N to 30N and from 9OW westward 





82 


to 855. These updated analyses are the ones which are saved in the historical 
files. About 500 radiometersondes are used in this final analysis. 

Based on the findings of this study, several modifications were made to 
the model which are: (1) the new calculations were confined to the troposphere, 
i.e., the layer from 100 to 30 cb; (2) for convenience and economy relative to 
computer costs, the geostrophic wind was used in place of the stream function 
wind; (43) based on Gates' (1960) tropospheric data for January and the results 
shown in Fig. 5, (the dashed line is to the left of the solid line in the tropo- 
sphere), the static stability parameter was set as a linear function of pres- 
sure. An artificial "wall" also was placed in the computation in that the 
values of the sine and Coriolis parameters at 4N were assigned to all grid 
points south of that latitude. 

Some preliminary results of this new study are shown in Tables 8 through 
14. The monthly average of the diabatic processes for the layer from 100 to 
Bomepeas a function of latitude is listed in Table 8. The transition from mean 
heating to mean cooling occurs near 50N. Note that the values and sign south 
of 4ON apparently support the criticisms of the NMC analyses relative to the 
boundary problems. Maximum heating in the layer occurs near HON while maximum 
cooling exists near 65N. Although Davis (19643) indicates that the volume be- 
tween 2ON and 7ON should undergo net cooling for January, this new calculation 
shows that for 1969, the net effects of latent heat release, solar insolation, 
and the transfer of sensible heat across the interface dominated the long wave 


ecoling south of 50N. 


Table 9 contains the monthly average of the generation of zonal APE as a 
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TABLE 8 


MONTHLY AVERAGES OF THE DIABATIC PROCESSES IN THE LAYER 
FROM 100 TO 30 cb AS A FUNCTION OF LATITUDE IN THE UNITS 
10° kj sec”* ton™* FOR JANUARY 1969 


Latitude 

5.0 2.4 
ie ere 
10.0 2a 
12.5 1 
15.0 0.5 
IS 0.5 
20.0 O52 
22.5 O-2 
25.0 Oe 
27.5 0.2 
Ole 0.6 
B25 eis 
35.0 0.9 
47 5 el 
Ie. 1 
42.5 eu 
ngs 0 10 
47.5 0.5 
50.0 -0.0 
52S -O.4 
55.0 On 
57.5 =e 
60.0 -1.1 
62.5 -1.2 
65.0 -1.4 
67 .5 alee 
FESO) -1.4 
72.5 =1L. 8 
75.0 -1.4 
77 .5 -1.0 
SOaC =0:.6 
82.5 -0.7 
85.0 -0.6 
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TABLE 9 


LATITUDINAL VARIATION OF THE GENERATION OF ZONAL 
AVAILABLE POTENTIAL ENERGY IN THE LAYER 100-30 cb 
COMPUTED FOR THE RING CENTERED AT THE INDICATED LATITUDE 
IN UNITS 107* kj m= sec™* FOR JANUARY 1969 


Latitude G pier eaee 
Z Deviation 

pire o. DP ee. Die 
7.5 4o5 .8 293.7 
TOO PE GPIASS 155.7 
12.5 100.2 96 .6 
15.0 Om 65.5 
ila 4.2 59.6 
One 22.0 54.8 
2225 = or 60.1 
25.0 7.6 38 .9 
ea o's - 5.6 40.8 
70). © Leen Beye lk 
S25 Ie AS one 
a 5-0 0.3 Sa 
47.5 -17.4 PPS 
hO.O -55 .8 50.4 
NAS oom 68 .2 
45.0 -{{ 3 Saal 
Home 5 “45.6 98.9 
50.0 Le ah es 5 
52.5 5 5 164.0 
Si 410 IES) 1S 0 
Sioa 197 .4 1g AL 
60.0 eye oh G2 AL 
62.5 ZO .6 7 a6 
65.0 287 .9 152 27 
645 euGae 18h .9 
70.0 25.2 eli. S 
EIS) PCS FEO 6 
75.0 BOO 192.0 
TENS: 220 206.0 
80.0 164.6 265.2 
82.5 13925 280.9 
a520 IyeGe, S20: / 
87 .5 15256 461.2 
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hoMmetion of latitude bands flor the layer from 100 Go 30 ch. Ss conpartseaswid 
Table 4 shows that the same pattern exists for the comparable bands although 
January has a wider transition zone and larger values of generation. This is 
not surprising since the sun's declination places it well south of the equator 
for January and this month probably represents wintertime extremes for both 
diabatic processes and temperature distribution in the Northern Hemisphere. 
The new information in Table 9 is the comparatively weak generation in the 
tropical latitudes south of 30ON which again may be due to the sun's position 
throughout the month. Results south of 1ON are probably unrealistic since 
there is no reason to expect that a quasi-geostrophic model produces valid com- 
putations at low latitudes with their correspondingly low values for the sine. 
Table 1O displays the monthly average of the generation of eddy APE as a 
Riebron Of latitude bands for the layer from 100 to 50 cb. Again, a compari- 
son with Table 3 shows a similar pattern with the notable exception that the 
region near 65N is no longer a maximum during January. The area associated 
with the polar frontal activity from about 35N to 50N is the maximum zone of 
destruction of eddy APE for this month. Latitude bands south of 30N do show 
that the rate of destruction of eddy APE decreases and changes to weak gener- 
ation south of 20N. (The area in the latitude band centered at 15N is larger 
by about 1.25 compared to the area in the latitude band centered at HON.) 
Table ll shows the layer variation of the generation of zonal and eddy 
APE for two domains. The larger region, from 11.25N to 88.75N, was chosen as 
being that region in which the quasi-geostrophic theory is reasonably appli- 


cable and the other region was selected so as to offer a comparison with March's 
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TABLE LO 


LATITUDINAL VARIATION OF THE GENERATION OF EDDY 
AVAILABLE POTENTIAL ENERGY IN THE LAYER 100-30 cb 
COMPUTED FOR THE RING CENTERED AT THE INDICATED 
LATITUDE IN UNITS 107* kj m~= sec7+ FOR JANUARY 1969 


Pelbaiscuele G eee 
e Deviation 

5.0 25.5 50.8 
WS) eal el .2 
1h © 20.8 1225 
1255 9.2 LS 
15.0 22.9 1S 
ly 35 ae 9.0 
2020 = 5.9 se, 2: 
Ly Ns = ed 15.7 
25.0 =a ge 
27 5 - 15.9 2S 
30.0 - 2).3 22.1 
32.5 - 41.7 Ale) iS) 
45.0 - 60.6 howe 
Die - 62.3 baie 
hO.O - 92.1 Teoma 
Nea =101.0 60.2 
h5.0 Sih oe 62.5 
“7.5 Sells ben 
BO.0 - 69.0 So 
Sel a's - 42.2 50) 
55.0 - 18.4 hh .g 
Pie - 6.9 39 .O 
60.0 jleezt alee 
62.5 - 4.6 Seg 
65.0 one BSc 
eS =) ioe Woe 
7 Ono = 0.6 290.6 
72.5 = 05 2Cnm 
75.0 = 7.6 See 
77.5 - 16.5 38.8 
80.0 - 16.8 Borne 
82.5 - 10.4 2S <5 
85.0 - 6.6 16.2 
Hs - 2.2 8.3 
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TABLE 11 
LAYER VARIATION OF THE GENERATION OF ZONAL AND EDDY AVAILABLE 
POTENTIAL ENERGY FOR THE INDICATED REGIONS IN THE UNITS 
107° kj m= sec7* cb~+ FOR JANUARY 1969 


Ie 5N = Gon oll 


Layer S standard G. Standard 
Deviation Deviation 
100 - 85 eb S10) 2 bie -54.0 alg 
85 = 70 cb 82 .3 an -53.5 18.4 
TO =- 50-eb 75 4 40.6 -42.5 19.5 
50 - hO ecb Fain 46.9 -2/.9 20.9 
4O - 430 cb 80.0 55.9 257-50 43.6 


20. (oN = 6875 


100 - 85 eb 105.9 Wh A -85.5 hO.O0 
65) =) 7 emeb 105 .6 30.9 -79 .8 28 .6 
TO wen Gi) ele Coyh 3 SS -h7 9 28.3 
50 - hO eb 95.4 hg. -O.) Bent 


4hO - 30 eb OT Oo 55.2 -h9.6 We 5 


results. A comparison with Table 4 indicates that aside from the intensity 
differences, the pattern is the same. For the smaller region, about 44% of 
the total contribution to the generation of zonal APE occurs in the lowest 40 
cb while slightly more than 70% occurs in the lowest 50 cb. Likewise, almost 
60% of the contribution to the destruction of eddy APE occurs in the lowest 40 
cb and almost 80% occurs in the lowest 50 cb. Values very close to these same 
percentages apply to the larger domain. Note that there is not a near balance 
between the zonal and eddy rates in the lowest layer for January. The results 
also suggest that within the tropics south of 30N, processes responsible for 
the generation of zonal and eddy APE are less intense at all levels during the 


wWillver . 
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Table l2 contains the contributions from the first fifteen waves to the 
generation of eddy APE for the two different volumes. In each domain, wave 1 
is the most destructive of eddy APE on the average while wave 4 is a secondary 
maximum. Shorter waves appear to be insignificant but the previous remarks 
concerning this scale still apply. The standard deviations show that wave 2 
displayed considerable variability for the month. Including the tropics dimin- 


ishes the net destruction of eddy APE at nearly all of the scales presented. 


TABLE le 


HARMONIC ANALYSIS OF THE GENERATION OF EDDY AVAILABLE 
POTENTIAL ENERGY IN THE LAYER 100-30 cb FOR THE INDICATED REGIONS 
IN THE UNITS 10°* kj m7* sec™* FOR JANUARY 1969 


eee WN SEG Sian s 26 Noo a, oll 
Average Standard Average Standard 

Number ane eres 
Deviation Deviation 

il ~6.2 5 cil -9.2 S20 

2 aS 6.9 25.1 10.4 

4 oe) SAE: -h .6 52 Ik 

h =5.5 hh -8.3 6.8 

5 -4.4 5c nS yenen 

6 “2.1 ey th “2.9 9.2 

i" =4.6 Bin ~-h 9 WS 

8 -1.4 ies -2.0 Bat 

9 =OnS 1.4 -0.8 2.0 

10 -O.1 Pal = 0929 | De 

iLL =-0.2 0.8 -0.3 ne 

12 =0519 O25 a0) 2 One, 

Ls sO I. 0.5 -0.2 ome. 

14 ove OS 0.0 0.6 

15 0.2 O.4 Cull 0.6 


Finally, Table 13 lists the monthly average of the generation of zonal and 
eddy APE for the different regimes. For the volume between 100 and 40 cb and 


between 11.25N and 88.75N during January 1969, the average generation of zonal 
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TABLE 13 


MONTHLY AVERAGE OF THE GENERATION OF ZONAL AND EDDY 
AVAILABLE POTENTIAL ENERGY FOR THE LAYER 100-30 cb 
BOUNDED BY THE INDICATED LATITUDE CIRCLES IN THE UNITS 
107* kj m7 sec”™* FOR JANUARY 1969 


11.25N - 88.75N 


G Standard G Standard 
‘ Deviation = Deviation 
54.9 19.6 19) 5 1 T2506 


28.75N - 88.75N 
ee 222 -44.4 eee 


APE was 55 x 10°* kj sec’* m= with a standard deviation of 20 units while eddy 
APE was destroyed at a rate of - 29 x Oza kj sec”* m7" with a standard devia- 


imonmot Ie units. 





CHAPTER 5 


SUMMARY AND SUGGESTIONS 


See CONCLUSIONS 

This study has examined results of the calculations of diabatic processes 
and the generation of available potential energy. Computations were made over 
the Northern Hemisphere between 18.75N and 88.75N in the layer from 100 to 10 
cb for March 1963. The input data consisted of the grid-point values from the 
NMC daily objective height analyses at eight levels within this layer. Dia- 
batic processes were calculated using a seven-layer quasi-nondivergent model in 
which the lower boundary processes and the vertical changes of the total deriva- 
tive (horizontal) of geostrophic potential vorticity are computed. As a result, 
fiomecie t1rst time, much better resolution of the vertical structure of ihe 
diabatic processes and the generation of APE, based on observational data, is 
available. The model further includes a more realistic treatment of the lower 
boundary effects. Specifically, the variable terrain height data due to 
Berkofsky and Bertoni (1955) and Cressman's (1960) drag coefficients are used 
to calculate the lower boundary vertical velocities and frictional effects. 
Also, in accordance with quasi-nondivergent theory, static stability is allowed 
to vary with pressure. 

Values of the calculated diabatic processes south of 30N indicate that 
the boundary problems associated with the analysis procedure using the NMC 
Octagon produced unreliable results. Also unusable were the computations above 


40 cb due to the failure of the model in the high troposphere. Problems arise 
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in this region because: (1) of the existence of a maximum wind layer which 
leads to a significant imbalance between the local time change and the hori- 
zontal advection terms involving geostrophic potential vorticity, and (2) the 
value of the static stability parameter becomes much larger here than in the 
lower troposphere. Vertical summation then propagates the failure throughout 
tme higher layers. Since the product of the square of the Rossby number with 
the Richardson number must be of order one in quasi-nondivergent theory, the 
existence of a jet or a very stable region invalidates this relationship and 
likewise the theory. Thus only the results in 18% of the earth's atmosphere 
are available for analysis. 

Monthly mean values of the diabatic processes in the layer from 100 to 40 
cb indicate that the change of surface cover to permanent snow or ice-cover is 
quite effective in separating the mean heating region from the mean cooling 
region in the troposphere. For March 1963, this separation occurred at about 
53N and extended almost vertically above this latitude throughout the tropo- 
sphere. The averages for the individual layers reveal the complexities involved 
when latent heat release, solar absorption, and the transfer of sensible heat 
from the surface to the atmosphere act to overcome the cooling effect of the 
long wave radiation component. In particular, although it would seem that the 
long wave component probably should dominate the layers above 50 cb everywhere, 
the heating patterns found at 40 cb suggest otherwise. Apparently the existence 
of regions of high cloudiness is required in order to trap the long wave emission 
and thus warm the layers beneath the cloud bases. 


Using the daily values of the diabatic processes in the four layers from 
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100 to 30 cb, the daily generation of APE was calculated. For March 1963, the 


- - -1 
average generation of zonal APE was 45 x 10 : kj m = sec. and eddy APE was 


destroyed at the rate of - 22 x 10 * kj m ae (These values apply to the 
troposphere which covers only one-half of the Northern Hemisphere's area.) 
Therefore, heating (cooling) takes place in lower (higher) latitudes and warm 
(cold) eddies are cooled (heated) everywhere on the average. Tropical lati- 
tudes and in fact one-half of the area of the Northern Hemisphere were not in- 
cluded because of the data limitations. It is probably reasonable to conclude 
that the tropics produce a strong generation of zonal APE. But due to the lack 
of significant wave structure similar to that found in the middle latitudes, 
little contribution to the eddy component should be expected on the resolvable 
Seale. Based on this reasoning, for the Northern Hemisphere, the generation of 
zonal APE is probably larger and the destruction of eddy APE is probably smaller 
than the values given here. 

In the region that was analyzed, on the average most of the generation of 
zonal APE and the destruction of eddy APE occurred in the lowest 40 cb of the 
troposphere. Further, in the lowest 15 cb there was a near balance between the 
zonal and eddy components. This suggests that most of the net production of 
APE takes place above the boundary layer. Although the net result in the 
troposphere was destruction of eddy APE, the layer from 50 to 40 cb did show a 
small average generation for the month. A strong average generation between 
4ON and 50N in this layer was responsible and this generation was probably due 
to the baroclinic activity along the polar front. High clouds in the warm sec- 


tor and strong cooling in the clear and cold regions would be very effective in 
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generating eddy APE. At the same time, the largest destruction took place 
throughout the troposphere in the baroclinic zone associated with the arctic 
front near 65N. Cooling of the warm eddies is probably the key mechanism there. 
Harmonic analysis reveals that the waves from 2 through 8 are responsible 
for the majority of the destruction of eddy APE. Very little energy transfor- 
mation occurred in the short waves. But since the observational data are not 
adequately distributed, the exact role of the short waves remains indistinct. 
Although wave 1 had little contribution in the troposphere on the average, 
strong periods of both generation and destruction were encountered throughout 
the month at this scale. It was particularly strong in helping to produce a 
net positive result in the layer from 50 to 430 cb. Waves 2 and 4 were the 
dominant modes within the troposphere and indicate that the scale of the land- 
sea distribution and the thermal contrasts were very significant for this month. 
The stated purpose of this investigation was to determine the feasibility 
of using a multi-layer quasi-nondivergent model to calculate the diabatic 
processes and the generation of APE. It has been demonstrated that the model 
is not applicable in the high troposphere and the stratosphere. However, compu- 
tations indicate that most of the processes which substantially affect the de- 
sired information occur within the troposphere itself. Based on these results, 
there seems to be little doubt that given adequate information relative to the 
mass-structure of the atmosphere, the model is more than satisfactory in de- 
tailing the vertical structure of the diabatic processes and the intensity of 


the general circulation within the troposphere. 
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5.2. SUGGESTIONS FOR FUTURE RESEARCH 

A very serious problem confronting observational studies is the lack of 
adequate data on a global basis. The problem is no less serious at the 
hemispheric level. But it will be some time before this situation improves. 

At least the available data should be handled in such a way so as to extract 
the maximum value. It is well within the present capability of numerical 
analysis to produce results which, although they may not be right everywhere, 
at least exhibit vertical consistency. In addition, the continued imposition 
of a lateral boundary well into the area of computational interest makes little 
sense at this date. In order to avoid the problems encountered by this study, 
the calculations should be made using analyses which cover a larger area and in 
fact all of the Northern Hemisphere. The analyses produced at FNWC are ideally 
suited for this task. Further, numerical analyses should have some way of 
showing the distribution and amount of input data. Perhaps, packed with the 
data on the magnetic tape, the number of reports used in the analysis would 
suffice. Then, knowing the standard distribution of the reports and the average 
number reporting, a more critical evaluation of the results could be made. And 
finally, it remains to be determined exactly what constitutes the "best" 
analysis of the mass-structure of the atmosphere. 

The use of the variable lower boundary values at a fixed pressure level 
was made in this study. Perhaps, as a result, the near balance in the zonal 
and eddy components for the generation of APE was calculated in this layer. It 
ls true that the number of grid points at which the pressure at terrain height 


is lower than 85 cb is small. Nevertheless, a comparison utilizing these dia- 
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batic effects as a function of terrain-height pressure should be made. 

Explaining the failure of the model in the high troposphere is based on 
showing that an imbalance exists among the terms. This imbalance occurs in two 
ways. The expression for calculating the diabatic processes above the lower 
boundary, (2.22), can be rewritten as 


p 


O 
= — + by p P 
HL [H oP, Bp |] Gere 


It was shown that a significant change occurs in the relationship of the terms 
in B, i.e., the local and advective changes in the geostrophic potential vor- 
ticity, in the high troposphere. This imbalance also affects the relationship 
between the processes at the lower boundary, Ho and the integrated effect of 
B. A serious imbalance at the higher levels of the model is made more so by 
the fact that the multiplicative factor, p/P» increases by an order of magni- 
tude there. Now obviously, if this imbalance could be made smaller, then the 
calculations might be more reasonable. How to do this is another problem. As 
can be seen in (2.23), the behavior of the static stability parameter does not 
affect the calculations of the generation of APE since H/( op) is involved there. 
But suppose that the imbalance of the given terms is real and that ne- 

glected processes act as the cancelling agents. Formulation of the model is 
based on the vorticity and thermodynamic equations. No approximations of any 
consequence are involved in the thermodynamic equation. The important assump- 
tions enter through the simplified vorticity equation. Since vertical veloci- 
ties tend to be smaller in a more stable regime, neglecting the vertical advec- 


tion of horizontal momentum would not appear to be serious near the jet stream. 
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Wiin-Nielsen (1959a) has shown that when horizontal advection is calculated using 
a nondivergent wind, absolute vorticity must be replaced by a constant value of 
the Coriolis parameter in the divergence term. This is done so as to avoid a 

net production of thermal vorticity. But in a diagnostic investigation, this 
point is irrelevant. Now it is known that the magnitude of the shear relative 
vorticity is large near the jet. In fact, regions of negative absolute vorticity 
have been computed on the anticyclonic side. Assuming that cw/ dp remains small 
so that the advecting wind can still be considered to be horizontally nondiver- 
gent, what would be the effect of allowing absolute vorticity to have its full 
range in the divergence term of (2.1)? From (2.8), it is clear that small 

values of absolute vorticity could cause computational problems, for example, 

a value of zero. It is not immediately clear what effect such a change might 
have. But a lower limit could be placed on the calculated absolute vorticity 
value and some computations might then be evaluated. However, the system would 
then be energetically inconsistent. <A consistent set of equations allowing the 
use of absolute vorticity would probably be more complicated. 

Perhaps another important and neglected process is frictional dissipation 
in the free atmosphere as proposed by Kung (1966, 1967), for example. Presum- 
ing that the sense of his argument is correct, then this might be quite impor- 
tant in the regions near the jet stream. But if it is, then the invalidity of 
the quasi-nondivergent theory is merely reemphasized. 

Tt very well may be that no simple theory exists by which the diabatic 
processes in the higher atmosphere can be determined diagnostically. But by the 


very nature of the system, the only significant component in these regions 
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should be long wave radiational cooling. The increasing data being made avail- 
able through the efforts of the aeronomers may eventually allow computations to 
be made in this regime. Although the present study implies that most of the 
important results occur in the troposphere, the inquiry should be expanded as 


soon aS practicable to include the stratosphere and to verify its exact role. 





APPENDIX 


FINITE DIFFERENCE FORMS OF THE EQUATIONS 


A.1. BASIC INFORMATION 

The finite difference forms of the equations used in the calculations for 
the diabatic processes are provided in this appendix. All standard symbols 
are explained in the list near the beginning of this paper. Figure 40 shows 
the computational grids required to locate the necessary points. The upper 
diagram is required only for the vertical differentiation needed in calculat- 
ing the geostrophic potential vorticity while the lower diagram applies to the 
horizontal operators. Also, the relevant equation number from the text is 
given where it is appropriate. 


As for the finite differencing scheme itself, consider some f(x) where 


Ax* 


- tiene (a) 


ti eeea) ee tC) t "((5c) Axe) 


Subtracting the two series represented by (1) leads to 


PC PAK) =) fC ee) : 


Te O(Ax)= , (2) 


= 


which is a centered difference formula for the first derivative. The order of 


=. is a measure of the accuracy of the approxi- 


the remainder, in this case (Ax) 
mation. Usually, the higher the order, the better the approximation is. 


Adding the two series in (1) gives the centered difference form for the 


second derivative, 
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Fig. 30. Vertical and horizontal finite differencing grids. 
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E(x) = Teg) a PaCS, + f(x - Ax) e O(Ax)= ; (4) 


Ie 
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putational instability associated with a time-stepping or prediction procedure 
is not a problem. But since the centered differences have remainders of higher 
order than either forward or backward differences, they are to be preferred and 
will be used. Round-off errors due to the physical storage capacity of the com- 
pucer and truncation errors due to the differencing scheme can be problems. 
However, nothing can be done about the former and so long as Ax is sufficiently 
eieli when compared to the characteristic horizontal scale of the problem, the 


latter is negligible. 


A.2. GEOSTROPHIC POTENTIAL VORTICITY, (2.6) 


lites seostrophie potential vorticity is calculaved using 


e 2 
=. a 42 ee ce Pat Bree, Oe eo, (4) 
re BPSD IND, tO Ne 
5 2 = ef + ae = 
Mee? = 0 tO, a * Og Py 


A.%3. LOCAL TIME DERIVATIVE OF GEOSTROPHIC POTENTIAL VORTICITY 
Using a centered finite difference over 24 hr, the local time change of 
the geostrophie potential vorticity is calculated from 


0b ee | Sie 


at 2~Sti“‘<i«éiC GK OF OS | (5) 
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Likewise, all other quantities in (2.8), although capable of being computed at 


observation time, will be averaged over a 24-hr period using 


O,,, 720 oU 
(), = i - 


t 4 


A.4. STREAM FUNCTION, (4.1) 


First, rescale (4.1) using 


ei 
V = — VY 6 
: (6) 
Then, 
-  VzZ 
aS es (7) 
ona 
. ae ae ao ees es "4d iad 
Vu = 2f z + . . (8) 
ee | eC Cement) ea 
Lely inet eg) abst i,jtl i,j ee 
A.5. HORIZONTAL ADVECTION 
Advection using the stream wind is computed using 
(V-9{}) = K ve (0 Cane. 
“V a =e a ta = 
hq@f Lt ee ae ee ase) a 
(9) 
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A.6. DEPARTURE VALUE AT THE TERRAIN HEIGHT 
Consider a grid point at which the terrain height is related to the heights 


of two bounding pressure levels such that 


7 Za, (10) 


Letting D represent the difference between the standard height and the actual 
Meme OL a pressure level, called the departure value, then the specific tem-— 


perature anomaly for the bounding layer can be written 


SC =3 a (a) 


Do. = DD = StZo 27 (12) 


Z = z -D : (ig) 


Therefore, (12) becomes 


DD. - ST(Z - 
U ( ie 4s) 


ee i) teres cee 
yoy 1+ST (14) 
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A.{7. PRESSURE AT TERRAIN HEIGHT 


From the results of (14) and knowing Z,. from (13), using an interpolating 
O 


polynomial, the pressure at the terrain height is calculated from 


_ Ly(etby)* + ule tau*)] 
P 1 (a5) 

where 

a = 724.64 

b = 4.5562 

C= eit 

jm 3.272 

7, “2On7 saloe 
Uu = Po 
1.5 110s 

and 

y= low 


A.8. TEMPERATURE AT TERRAIN HEIGHT 


By definition, the specific temperature anomaly is 





Since the standard temperature at a given pressure level is calculated from 
T = 288. - 6.5 x10°%-2 ; (17) 


the temperature at terrain height is computed from 
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ae (288. - 6.5.x 10-2 eee eR (18) 
O 


A.9. GHEOSTROPHIC SURFACE WIND 


The geostrophic surface wind can be represented by 
> > 
VS hse (a) (19) 


Now the gradient of D at the terrain-height pressure requires the use of (12) 


but in this case a is known. Letting the subscript op represent these values 
O 


at terrain-height pressure, the geostrophic surface wind is determined from 


a5 
2 


V = ng_ 2 


D -D f 
O Pia L Ose les OY hk A 
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A.10. GEOSTROPHIC SURFACE VORTICITY 

From the same considerations used to obtain (20), the geostrophic surface 


vorticity can be computed from 


Peele FRICTION EFFECTS, (2.15) 
Combining (15), (18), (19), and (20), the frictional effects are computed 


from 
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where the brackets contain the terms inside the brackets in (20). 


A.12. OMEGA IN THE LOWER BOUNDARY, (2.19), (2.20) 
The omega due to terrain is calculated using a form similar to (9) except 
that the geostrophic surface wind is used in place of the stream wind. Cal- 


culating the omega due to friction is just a matter of multiplying (22) by 


-Ap/f. 


A.13. HEATING AT THE LOWER BOUNDARY, (2.17) 

Computing the lower boundary heating is straightforward using (15), (18), 
(19), and the vertical velocities. All the terms are averaged over the 24-hr 
period since the local time change term is also involved. The averaging is 


done using the form given in Section A.3. 
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